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Abstract 

Upon arrival to a ticket queue, a customer is offered a slip of paper with a number on it - indicating 
the order of arrival to the system - and is told the number of the customer currently in service. The 
arriving customer then chooses whether to take the slip or balk, a decision based on the perceived queue 
length and associated waiting time. Even after taking a ticket, a customer may abandon the queue, an 
event that will be unobservable until the abandoning customer would have begun service. In contrast, a 
standard queue has a physical waiting area so that abandonment is apparent immediately when it takes 
place and balking is based on the actual queue length at the time of arrival. 

We prove heavy traffic limit theorems for the generalized ticket and standard queueing processes, 
discovering that the processes converge together to the same limit, a regulated Ornstein-Uhlenbeck 
(ROU) process. One conclusion is that for a highly utilized service system with a relatively patient 
customer population, the ticket and standard queue performances are asymptotically indistinguishable 
on the scale typically uncovered under heavy traffic approaches. Next, we heuristically estimate several 
performance metrics of the ticket queue, some of which are of a sensitivity typically undetectable under 
diffusion scaling. The estimates are tested using simulation and are shown to be quite accurate under a 
general collection of parameter settings. 


1 Introduction 

In many service settings, newly arriving customers are given information about the number of individuals 
preceding them in line, even when this line is virtual. There are several ways in which this information may 
be passed on. For example, a customer visiting either a delicatessen or the department of motorized vehicles 
(DMV) is offered a ticket with a number on it and, via some physical display, is informed of the current 
customer being serviced. In restaurants, dinner parties may either be told about the estimated wait or told 
how many similarly configured dinner parties are ahead of them. Often these two forms of information are 
roughly interchangeable: given the service rate, knowledge of the queue length yields an estimate of the 
delay, and vice versa. Being informed about the delay in service provision, customers then choose whether 
to join the queue or to balk. 

The fact that a customer initially accepts the estimated delay and joins the queue does not guarantee 
that the customer will wait around until service can begin. Customers may renege on their initial decision 
and abandon the queue. In environments where customers are physically waiting in line for their service - 
such as at a bank or grocery store - abandonment is immediately apparent to service providers and other 
customers alike. However, neither the delicatessen service personnel, fellow ticket holders, nor potential 
ticket holders are aware when someone has chosen to abandon their ticket. This event is not discovered until 
the ticket’s number is called and no one responds. In general, the number of outstanding tickets may be 
larger than the number of customers actually waiting for service. 

A ticket queue refers to the setup typically employed in delicatessens, but can be thought of more generally 
as a mechanism for tracking the number of potential customers yet to be served and for maintenance of a 
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Figure 1: Queueing dynamics of the ticket and standard queues when abandonment occurs. 


first-come-first-served protocol. By potential it is meant that these customers have been triaged, joined 
the queue, and have committed in principle to be served, yet may ultimately abandon before service can 
actually begin. Some of the customers may renege on this implicit commitment and leave; reneging customers 
typically will not inform the system manager of their decision to forgo their place in line. As a result, what 
is generally perceived as the queue length the number of potential customers to be served will in fact be 
an upper bound. 

A visual comparison of abandonment in the two queueing types is depicted in Figure 1. The left side of 
the figure shows a progression of three events as they are experienced in the standard queue. The first image 
has a queue with four customers, numbered 1 through 4 and with customer #1 currently in service. In the 
next frame, customer ^3 abandons the queue, the system immediately detects the defection and customer 
#4 replaces customer #3 in line. Then customer #1 completes service. Finally customer $4 reaches the 
front of the queue and begins service in the 4th frame after the service completion for customer #2. 

The right side of the image depicts the same sequence of events but as they are experienced in a ticket 
queue. Notice that when customer #3 abandons his ticket, the system is unaware. Thus, the perceived 
queue length (including the customer in service) is four. After customer #T completes service, all customers 
advance and customer #2 begins service. After customer #2’s service is complete, an attempt is made to 
handle customer ^3, and it is determined immediately that customer ^3 has left. So customer $=4 begins 
service at this time. What we see here is that in two of the four frames, the ticket queue’s perceived length 
is larger than that of the standard queue. 

Several questions emerge: What are the dynamics of the collection of outstanding tickets? What is the 
actual number of customers remaining, i.e., those who have yet to renege? What fraction of customers 
renege? What fraction of customers balk? Does the difference in implementation of ticket and standard 
queues lead to a marked difference in performance between the two customer organizing methods? 

Some heuristics seeking to address these questions were provided by [? ]. Our paper revisits the ticket 
queue to provide new heuristics and intuition that is complementary to that of [? ]. In particular, informed by 
a heavy traffic limit theorem, we conclude that for highly utilized systems with relatively patient customers, 
ticket queues and conventional queues are indistinguishable. One of the conclusions of [? ] is that for 
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heavily loaded systems with relatively impatient customers, the ticket queue experiences a higher percentage 
of balking. In addition to these divergent insights, what sets this paper apart from [? ] is that it provides 
estimates of the distribution of the queue length process, which can then be used to estimate other stochastic 
elements. More details about these estimates are provided below. Moreover, general assumptions about 
the abandonment distribution and balking distributions are employed; random variables in [? ] are all 
exponentially distributed. 

The paper appeals to heavy traffic limit theory in the pursuit of its approximations and comparisons of 
the ticket and standard queues. In this regard, it is similar in approach to [? ], which studies a generalized 
diffusion-scaled single-server queue with abandonment or balking. Our paper concerns a realistic situation 
where there is both balking and abandonment. Moreover, technical complications arise as we base balking on 
the queue that customers perceive, whereas abandonment is based on the delays that customers experience. 
Ultimately, what one obtains under diffusion scaling is a regulated (at zero) Ornstein-Uhlenbeck (ROU) 
process whose constant drift is related to the difference between the arrival and service rates and whose 
restorative drift involves the derivatives of the abandonment and balking distributions, both evaluated at 
zero. The diffusion limit of the critically loaded ticket queue is identical, despite the fact that for the 
ticket queue, abandonment of a ticket is not matched with a shortening of the queue until the ticket of the 
abandoning customer reaches the front of the queue. Further, given the same primitive random variable 
elements, the diffusion-scaled ticket and standard queueing processes converge together to the same diffusion 
limit, a stronger notion of process similarity. In other words, not only are the dynamics of the ticket and 
standard queues asymptotically similar, the processes are asymptotically coupled. 

As mentioned above, we estimate the distribution of the queueing processes employing heuristic inter¬ 
pretations of the limit theorem. Additionally, we estimate performance metrics such as abandonment and 
balking probabilities and the expected number of abandoned tickets in circulation. Each of these additional 
quantities is typically lost under conventional heavy traffic limit theory, but are manifest on a smaller - 
i.e., more sensitive - scaling. The approach for estimating the expected number of abandoned tickets in 
circulation provides the additional insight needed to compare the subtle differences between the standard 
and ticket queues. 

We should also mention that the differences between standard and ticket queues also highlights the 
different ways that delay information can be communicated to the customer. This communication between 
the service and its customers is important because customers will make their decision to wait or to leave 
the queue based on the information that they receive from the manager of the system. For example, see the 
following papers on research on queues with delay announcements and estimation [???????]. 

As for the assumptions used throughout, the paper also fits within the growing literature on queues with 
generally distributed abandonment distributions, specifically those that are not exponentially distributed. 
Some recent examples of such papers include [????????]. The first two of these are in a multiserver 
setting, whereas the others are in the single server regime. The last two papers use measure-valued processes 
to model system dynamics. 

The remainder of the paper proceeds as follows. Both the ticket queueing model and standard queueing 
model (with abandonment and balking) are presented in the next section. The main result, Theorem 3.1, 
is presented in Section 3. Section 4 contains the heavy traffic-inspired approximations for performance 
measures and some interpretations. The proofs of the main results are provided in Section 5. Extensive 
numerical results are presented in Section 6. Concluding remarks and extensions follow. 

1.1 Notation 

We conclude this introduction with notational conventions. Let 7Z denote the set of all reals, 7Z+ the set of 
nonnegative reals, and J\f the set of strictly positive integers. For a Polish space S, let D(7Z + ,S) denote the 
space of right continuous with left limits functions from 1Z + into S. The Polish spaces we consider here are 
1Z + and 1Z 2 + . 

For ease of exposition, quantities that are related to the ticket queue will be appended with the subscript 
T and those associated with the standard queue will have the subscript S . We use the subscript a as a place 
holder for either ticket (a = T) or standard (a = S ) queues. The standard hazard rate function h : 1Z K > TZ+ 
is the ratio of the density of the standard normal distribution to the tail of the standard normal distribution: 
h(x) = (f){x)/{l — $(£)) f° r every x £ 1Z. For enumeration purposes, we use the letters i,j , and k to represent 
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nonnegative integers. The letters q, r, s, and t are used to represent time. Typically, the symbols 6, e, rj are 
used to represent small positive real numbers. In contrast, the letters I\ and L are used to represent large 
quantities, predominately as upper bounds. 


2 The model basics 

In this section, we provide the primitive random variables for the queueing processes, describe the construc¬ 
tion of the ticket and standard queues, and discuss the intricacies of simulating the ticket queues. 

2.1 Random variables 

We construct the ticket and standard queueing processes using the same collection of random variables. 
Each customer that arrives after time zero has an interarrival time, (potential) service time, initial time 
tolerance and (potential) abandonment time. For the i th customer, these times are captured in the quadruple 
(u,i,Vi,bi,di). The letters ’b’ and ’d’ denote balking and deadline, respectively. The mutually independent 
sequences {ui,i > 1}, {vi,i > 1}, {h,i > 1}, {d t . i > 1}, which are all defined on the same probability space 
are each i.i.d. The exogenous arrival rate of jobs is A and the service rate is /i. The sequences 
{ui,i > 1} and {vi,i > 1} have unitary means and the interarrival and service times of the i t h customer are 
Ui/X and Vif/j,, respectively. The arrival time of the ith job occurs at time 

i 

u = (v a )X! u t 

3 =1 

The unitary interarrival and service times have variances a% and of, respectively. The quantities 6,; and di 
represent the random variables associated with balking and reneging, respectively. Let F), and F ( j denote their 
respective cumulative distribution functions. We assume the these functions both vanish at zero. Moreover, 
their derivatives exist at zero and the sum of these derivatives is strictly positive. We define the sum as 

0 = E'(O) + F'(O). (2.1) 


2.2 Initial conditions 

Consider jobs that are present at time zero. These jobs do not require arrival times. Neither do they require 
balking random variables. Hence, for these initial jobs, we provide only two sequences of random variables: 
unitized potential service times, {vi,i > 1}, and residual deadline quantities, {di,i > 1}. Let <3(0) denote 
the number of initial jobs. If initial job i < Q(0) has not begun service before time di, this job will abandon. 
If this initial job has not abandoned, then its service time will be Vi/fi. The initial potential service times 
are i.i.d. and have the same distribution as the potential service times of the jobs arriving after time zero. 
The residual deadlines do not necessarily have the same distribution. Let F,]j be the cumulative distribution 
of di . We impose the following uniform restriction on their distributions near zero: There exists an / > 0 
and an ho > 0 such that 

FdAh) , i w, , , 

sup — - < /, vh < hQ. 

i h 

The workload at time 0, denoted IF(0), is the amount of effort required to process jobs present at time 
zero. Because of deadlines running out before service has begun, some of the Q( 0) jobs present at time 
zero will not be served. Therefore, it is not as simple as adding up the service times of the first Q(0) jobs. 
Instead, let uq denote the cumulative amount of server effort required among the first i jobs in-queue at time 
zero. The i th job will be served if and only if it is sufficiently patient, or if di > u>i- 1 , where w o = 0. We 
can define the wf s recursively: 

i * ~ 

■Wj = ■ 1 (dj > U)j- 1 ) = Wi -1 + -JJ ■ 1 (di > 1 ), i > 1 . 

It follows that IF(0) = u>q ( 0 ). 
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2.3 The queueing processes 

Let the process Qt = {Qr(t),t > 0} track the dynamics of the ticket queue and Qs = {Qs(t),t > 0} track 
that of the standard queue. The difference between the ticket and standard queue is the timing of when the 
abandonment is accounted for. Otherwise the system dynamics are identical. For example, customers arrive 
to the ticket and standard queue in the same sequence, with the same balking and reneging tendencies, and 
with the same service requirements. The common arrival process is A = {A{t),t > 0}. 

All other processes are indexed by the abandonment protocol a £ {S, T}. In particular, B a = {B a (t),t > 
0} is the balking process that tracks as a function of time the number of arriving customers who leave 
immediately upon arrival. The reneging process R a = {R a (t),t > 0} tracks the number of jobs that arrive 
after time zero , who have abandoned, and whose abandonment has been detected in the system. Recall 
that in the standard queue an abandoned customer is immediately detected once they leave, whereas for 
the ticket queue abandonment is only apparent at the moment at which service for that customer would 
have begun. The process S a = {S a (t),t > 0} tracks the number of service completions, of jobs that arrived 
after time zero , as a function of how much effort the server has expended. Relatedly, the busy time process 
T a = {T a (t), t > 0} reports how much time has been spent processing jobs as a function of time t, including 
initial jobs. The idle time process I a = {I a (t),t > 0} is complementary to T a : I a (t) = t — T a (t) for each 
t > 0. Lastly, the workload process W a = {W a (t),t > 0} reports as function of time the amount of effort 
required by the server to process those customers currently in-queue that will not abandon. 

The system must clear out all initial jobs in-queue before it can start processing jobs that arrive after 
time zero. Let Q a = {Q a (t ), t > 0} track the number of remaining initial jobs at time t. Some of these jobs 
may abandon. Let R a = {R a (t),t > 0} track, as a function of time, the number of jobs who arrive before 
time zero, who have abandoned, and whose abandonment has been detected in the system. The standard 
and ticket queues start with the same collection of jobs at time zero. Hence, among these initial jobs, the 
same subset of jobs abandon both the standard and ticket queues. What is different about the abandonment 
processes is the timing of when the abandonment is detected. What is identical between the ticket and 
standard queues is the initial job service completion process: S = {S(t),t > 0}, which is defined in the next 
section. 

Equations governing the ticket and standard queue are below. For each t > 0 and a € {S,T}, 


Qa(t) = Qa(t) + A{t) - B a (t) - R a (t ) - S a ((T a (t) - W(0))+), 

(2.2) 

Qa(t) = Q(o) - R a (t) - S(t), 

(2.3) 

Aft) = sup | j > 0 : ^ m /A < 

(2.4) 

A(t) 

Bait) = ^2 1(6* < Q a {ti-)/n), 

i= 1 

(2.5) 

Mt) 

Rr(t) = ^2 1 (h > Qt{U — )/a0 • 1 (d; < Wt(U—)) ■ 1 {WriU—) <t — tf), 

i— 1 

(2.6) 

Q(0) 

Rr(t) = ^2 < Wi- 1) ' < t), 

i—l 

(2.7) 

m 

Rs(t) = *^2 1 ^ i > Qs{ti-)/n) ■ 1 (di < mm(W s (ti-),t - U)), 

i—l 

(2.8) 

Q( o) 

Rs(t) = ^2 l idi ^ min(uij_i, t)), 

i= 1 

(2.9) 

Ta(t)= I l(Q a (s) > 0)ds, 

Jo 

(2.10) 
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and 


I a {t) = t- T a (t), 


(2.H) 


m 

w a (t) = W( 0 ) - T a (t) + £>*//*) ■ > Qa(ti-)/ij) ■ 1 (di > W a (u-)). ( 2 . 12 ) 

i -1 

Interpreting (2.2) and (2.3), the queue length process consists of the remaining initial jobs and may 
increase with each arrival, provided the corresponding customer does not balk. In addition to the departure of 
the initial jobs, the queue length process decreases whenever there is an abandonment or a service completion 
among the jobs that arrive after time zero. Initially, service is allocated entirely to the initial jobs and remains 
so until, W{ 0), when those jobs have departed entirely. At this point service allocation is given entirely to 
jobs arriving after time zero. The initial jobs experience only abandonment and service completion. The 
remaining job process is a decreasing function, hits zero, and remains there. 

As for the balking process and (2.5), an arriving customer joins if their initial delay threshold hi is 
sufficiently large. When a customer arrives to the system, it is first triaged and given the opportunity to join 
the queue. Under the ticket queue implementation, joining the queue involves the acceptance of the offered 
numbered ticket. In addition to knowing which ticket it will be given, the customer is told the number of the 
ticket holder currently in service. For the standard queue, joining the queue involves standing in a physical 
line. In both cases, the decision of whether or not to join the queue is based on the customer’s expectation 
of the delay until service and her tolerance for such a delay. Customers convert the queue information into 
an expectation of delay until service. We assume the conversion is naive and the same for both ticket and 
standard queueing environments: Given the queue length, the customer estimates the delay by dividing the 
queue length by the service rate /z, a quantity assumed to be known by all customers. That is, if the i t h 
customer arrives at time t, the customer joins the queue if bi > Q a (t—)/n ; otherwise the customer balks. 
The sequence of balking tolerances has common distribution function Ff,. Naturally, the probability that a 
customer arriving at time t balks is Fb(Q a (t—)/n). 

Consider (2.6) and (2.8). A customer who joins the queue is not guaranteed to stick around for service. 
If the delay that customer i experiences in the queue reaches di then that customer will abandon from the 
queue. The time that a customer arriving at time t would have to wait is captured by W a (t—). An analogous 
workload process is the main object of study in [? ]. The distribution of the abandonment time random 
variables di is denoted F ( f, the ‘d’ stands for customer deadline. It follows that a non-balking customer 
arriving at time t will abandon the queue with probability Fd(W a (t— )). What sets the ticket queue apart 
from the standard one is the time at which the process Q a reflects the abandonment of a job. Hence the need 
to express Rr(t) and Rs{t) separately in (2.6) and (2.8), respectively. Consider (2.6) and how it captures 
customer abandonment. Not only must the tolerance di be smaller than the delay that the customer must 
endure before service, the system does not know that the customer has abandoned until that delay has 
expired. Equations (2.7) and (2.9) are the analogous formulations for reneging customers who are present 
at time zero. 

The workload process is also referred to as the virtual waiting time process because it tracks, as a function 
of time, the amount of time a sufficiently patient, non-balking customer would have to wait before receiving 
service. The virtual waiting time process increases by the service time whenever a job arrives to the system 
that will eventually receive service. The process decreases at rate one whenever it is greater than zero, or 
equivalently, whenever the queue length is nonzero. The cumulative amount that the workload has decreased 
by time t is precisely equal to the total busy time T a (t). 

2.4 State space descriptors and simulation of the ticket queue 

Simulating the ticket queue is more complicated than simulating the standard queue. We describe below the 
intricacies of simulating the ticket queue under both Markovian and non-Markovian assumptions. We use 
simulation later to assess the accuracy of our approximations and heuristics. 

Under Markovian assumptions, the state space of the ticket queue can be captured by a vector of zeros 
and ones. The length of the vector corresponds to the number of outstanding tickets, including the ticket of 
the customer currently in service. (For convenience, assume that the first element of the vector is the leftmost 
element.) If the vector has a nonzero length, the first element corresponds with the customer in service and 
by convention is a one. The other elements correspond with the other unresolved tickets. Further, the order 
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of these elements reflects the relative order of the corresponding customers’ arrivals and, under our first- 
come-first-served assumptions, the order in which resolution will take place. Ones in the vector represent 
customers who have not abandoned the queue. Some of these customers may abandon before resolution 
takes place. When a customer abandons, the corresponding element turns into a zero. When service of a 
sufficiently patient customer takes place, the state vector shifts to the left because at least the first element 
of the vector must be removed. Either the second element is a zero or it is a one. In the latter case, the 
entire state vector shifts by one element. If the second element is a zero, this represents an abandoned ticket. 
Starting with this zero in the second element, all contiguous zeros will be will be removed from the state 
descriptor. The assumption here is that resolution of abandoned tickets is instantaneous. When a job arrives 
to the system and chooses to join the queue a one is appended to the end of the state vector. If the customer 
balks, no change in the state takes place. 

The system transition is governed by exponential clocks for each unresolved ticket that has yet to be 
abandoned and is not being processed, one clock for the job in service, and one for the next arriving job. If 
the clock associated with an unresolved ticket not in service expires, then the associated customer abandons 
and the element in the state descriptor changes from a 1 to a 0. If there is an arrival, then another random 
variable is generated and compared to the weighted queue length to determine whether the customer balks; 
if balking occurs the state does not change. If the clock associated with the job in service expires, service 
completion ensues and the state changes as described above. 

Alternatively, one could have one exponential clock for all unresolved tickets that are not in service 
and have not been abandoned. The rate of this clock is equal to the number of such tickets multiplied 
by the abandonment rate of a single individual. If this clock is the one that expires next then the actual 
abandoning customer is found by randomly choosing between the non-abandoned waiting customers with 
equal probability. When there is a change in system state, this clock must be recalculated because the 
number of unresolved tickets will have changed as well. 

Under general assumptions on the random variables, the state space must contain the residual interarrival 
time of the next customer to arrive, the residual service time of the customer at the front of the queue, and 
for each customer yet to reach the front of the queue, the residual abandonment time. There is an alternative 
state space formulation, if one is content with only knowing which of the customers in-queue will eventually 
be served. For this alternative, one must track the virtual waiting time process, Wt, which yields as a 
function of time the amount of time that a customer must wait until service begins, and the eventual service 
time of jobs that will be served. These service times can be kept in a vector, similar to the vector of zeros 
and ones above. From the time of their arrival, jobs that will have abandoned before reaching the front of 
the queue have a zero in their corresponding element of the vector. The virtual waiting time process jumps 
at the time of an arrival by the service time of the corresponding customer only if this customer actually 
joins the queue and is sufficiently patient; see for instance [? ] or [? ]. What is lost in this formulation is 
the timing of the individual jobs’ abandonment times. Gained is the freedom from having to track residual 
abandonment times for each job in-queue. 


3 Comparing the queue processes 

In this section we provide the main result, a heavy traffic limit theorem that serves as the theoretical 
underpinning of the heuristics forwarded in the subsequent section. 

To facilitate comparing the ticket and standard queue, we appeal to heavy traffic limit theory. To this 
end, we consider a sequence of systems, indexed by n. The arrival and service rates of the nth system are 
A" and /i n . Equations (2.2)—(2.10) have straightforward analogs with the A and ^ replaced by X n and /r n , 
respectively. For each t > 0, and a G {S', T}, 

Ql(t) = Q n a (t) + A n (t) B2(t) R n a {t) - S™(T™ (t) - W n m, (3.13) 

QZ(t) = Q n (0)-R2(t)-S n (t), (3.14) 

A n (t) = sup I j > 0 : ^ Ui/X n < 1 1 (3.15) 
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(3.16) 


A"(t) 

B n a {t) = E < Q n a W-)/n n ), a G {5,T}, 

*=1 

^(*) 

7^(t) = E 1 ^ > Qt(U-)/^) ■ 1 (di < W2(U-)) • 1 (W?(t<-) < t - ti), (3.17) 

i=l 

Q n ( 0 ) 

£?(t) = E ^ ^ ™"-i) • W-i < *), (3-18) 

A"(t) 

i2§(t) = E ^ > Qs(ti-)/n n ) • 1 (di < min(M^(ti—),i-*<)), (3-19) 

i=l 

Q n { 0 ) 

B n s{t)= E l(rf» < min^lLi.f)), 

Z=1 

W)= f im(s) > 0)ds, a G {S', T}, 

Jo 

A n (t) 

WS(t) = W n ( 0) - !£(*) + E («i/M n ) • l(6i > QS(*»-)/m") • l(d< > a € {5,T}, 

i=l 

and 

IZ(t) = t-TZ(t), ae{S,T}, 

where 

i 

*? = ( l/A")E«i. 

1=1 

i ~ 

W n (Q) =Wq H0) , m>o=°. and «>? = E WT ‘ 1 (4i > w"-i)> *>!• 

l=i M 

The associated scaled processes are Q” = > 0}, Q” = {Q”(t),t > 0}, A” = {A ra (t),f > 0}, 

B n a = {B2(t),t > 0}, i?2 = {J2S(t),t > 0}, = {^S(t),t > 0}, SI = {52(t),t > 0}, = {52(1),i > 0}, 

TZ = {TZ(t),t> 0}, 72 = 0}, and W" = {W^(t),t> 0}. 

We envision the arrival and service rates each being order n and differing by a quantity that is order 
\fn. So, in the absence of abandonment and balking, one would expect the queue length to be order y/n 
and for the workload process to be order 1 /y/n. In fact, this intuition is true in the presence of both 
balking and abandonment. Hence, we define for each a G {S, T}, the diffusion scaled queue length process 
Qa = {QaiOit — 9} and the inflated workload process W™ = {W£(t),t > 0}, where for each t> 0, 

Qa(t) = ^0- and W2{t) = y/nWZit). 

yn 

Notice that we do not scale time as the arrival rates and service rates are already proportional to n. Also 
notice that the balking and abandonment times do not change with n. The reasoning is that demand may 
change and service speed must adjust accordingly, however, individuals will still have the same desires for 
and assessment of service quality. 

We introduce the processes e2 = {e”(t), i > 0} for each a G {5, T}, where for each t> 0, 

/ A n (t) \ 

e n a (t) = K(t) + I ~ E < Q n a (TP-)/n n ) 1 • (3.24) 

The idea is to replace the reneging process 7?2 with a process that ignores whether the job has balked when 
considering whether it will renege. In reality, a balking customer leaves and, as a result, the question of 
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(3.20) 

(3.21) 

(3.22) 

(3.23) 



whether balking customers would renege is a moot one. The introduction of this process also eliminates 
the concern of when the reneging customer causes a decrease in the queue length. Here we assume that 
the customer never actually enters the queue. One further subtly is that the workload is replaced by the 
weighted queue length the quantity used to determine balking - so that one need only track the process Q ™ 
rather than the joint process (Q™, W™). Lastly, a benefit of this formulation is that it allows for the ticket 
and standard queues to be handled simultaneously. Ultimately, we will show that the process e™ is negligible 
under diffusion scaling, which partially argues why the diffusion-scaled ticket and standard queues converge 
together to the same limit. The process e™(-) also eliminates the reneging of the initial jobs altogether; the 
reneging of initial jobs is shown to be neglibile in Proposition 5.13. 

The process S™ tracks the number of customers (who arrive after time 0) served as a function of total 
effort dedicated to customers. As not all customers receive service, the service times that determine S™ are 
a subset of {vt,i > 1}, the collection of potential service times of arriving jobs. This subset differs for the 
ticket queue and the standard queue. The index of the i t h job whose service time contributes to S ™ - that 
is, who is actually served - is 


JaW = inf 



1 : ]T l(h > Q n a (t?-)/n n ) ■ l(di > WZW-) > i 


i > 1, ae{S,T}. 


The sequence of service times that are actually used is denoted {u”(i),i > 1}, where v™(i) = Because 

the service time of job j£(i) is independent of all random variables that dictate whether this service time is 
used, the filtered sequence > 1} is i.i.d. and has the same distribution as the original (unfiltered) 

collection of service times. Therefore, any property of the unfiltered service times - such as weak laws of 
large numbers or invariance principles - holds for the filtered sequence. Finally, we can write S™( t) for each 
t > 0 as 


5”(t) = sup 


k> 0:(l/M n )V>2(i) 



a g {. S,Tj. 


(3.25) 


The processes S = {S(t),t > 0} and S n = {S n (t),t > 0} are also defined analogously. 
Now we can write the diffusion-scaled queue length processes for each t > 0 as 


Q n a {t) = Q n ( 0) + A n (t) - M b (A n (t)) - M2JA n {t)) - %(t) %(t) - S2((T2(t) - W"(0)) + ) 


-9 / Q n a {s)ds + 


(A" - //>) 


t + Y2(t), ae{S : T }, 


(3.26) 


where, for each a £ {S,T}, Q™(0) = ( l/y/n)Q n (0) is the scaled initial queue length and, for each t > 0, 


M2,aV) 


A n (t) = (l/Vn) (A n (t)-X n t)., 
A n (t) = (1 /n)A n (t), 

[nt] 


(1/Vn) (nh < - F b (Vn.Q n a (t?-)/v n )) , 


2 = 1 
[nt\ 


(1 /Vn) E (l (di < Q n a (i?-)/n n ) - F d (V^Q n a W-)/n n )) , 


2=1 


S2(t) = (l/Vn) (S2(t)-n n t), 
e n a = (1 /Vn)C{t), 


(3.27) 

(3.28) 

(3.29) 

(3.30) 

(3.31) 

(3.32) 


%(t) 


i / \ r* 

~^= (F b (VnQ n a (t?-)/v n ) + F d (VnQ n a (t?-)/fi n )) - 9 ^ Q n a (s)ds (3.33) 

+ ^=(5”(t)- At ”min(t,lU"(0))), 
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and 


(3.34) 


K(t) = (£) W) = ((£) m. 

We refer to A n = {A n (t),t > 0} as the diffusion-scaled arrival process and to A n = {A n (t),t > 0} as its 
fluid-scaled analog. The reader may notice that the process 6™ = {d„(t). t > 0} has what looks like an 
instantaneous drift that is proportionate to the value of the scaled queue length process. The remaining 
processes are centered and diffusion scaled versions of their original analogs: Mff a = > 0}, 

MS,a = {M% a (t),t> 0}, 52 = {S2(i),i > 0}, cS = {<£(*),*> 0}, and Y£ = {Yff(t)’t > 0}. 


3.1 A heavy traffic limit theorem 

In order to prove a heavy traffic limit theorem, we assume for our sequence of systems indexed by n, that 
arrival and service rates are order n quantities and are asymptotically identical; that is, asn-> oo, 

\ n /n —y p, and p n /n —> fi. (3.35) 

Further, the difference between the two should be an order y/n quantity such that as we take the limit 
n —> oo, 

(A ra — n n )/\/n = (3 n —► ft £ (— 00 , 00 ). (3.36) 

One can refer to (3.35) as the heavy traffic condition; the expression implies that 

p n = \ n /p n _► 1, (3.37) 


as n —> 00 . We assume that the random variables associated with balking and abandonment are unaffected 
by the change in the index n. Define 

= h\] al + ct 2 s (3.38) 

as the standard deviation associated with the arrival and service times. Lastly, define B = {B(t),t > 0} as 
a Brownian motion with no drift and an infinitesimal variance of 1. 

The framework developed in [? ] justifies the alternative representation of (3.26): 

(QZ,Y:) = (<b e ,'M(Q n (0) + X2), (3.39) 

where (4>g, T#) : T>(TZ + ,TZ) D(IZ + ,'R,\) is a Lipshitz continuous map, X™ = {X™(t),t > 0}, and for each 
t> 0 and a £ {5, T}, 


X-(t) = A n {t) M£ a (A n (t)) - M£ a (A"(t)) - £(i) - 6™(t) - S2(72(t)) 


(a n _ zj n. 


(3.40) 


The elements of X n are those that either will converge to Brownian motions or that are asymptotically 
negligible. The limiting stochastic process is the following 

X = /3e + aB. (3.41) 


We now present our main result for the diffusion scaled queue length and workload processes. 

Theorem 3.1. If 

(Q n (0),W n (0)) =► (Q 0 ,Qo/fi), as n ^ 00 , (3.42) 

then 

((Q n s ,WZ,Yg)AQZ,W2,Y?))^((Q,Q/p,Y),(Q,Q/»,Y)), as n ^ 00 , (3.43) 

where Q(0) is equal in distribution to Q 0 , Q = 4>g(Q(0) + X), Y = 4'g(Q(0) + X), and together Q and Y 
obey the following stochastic differential equation 

dQ(t) = —9{f3/0 — Q(t))dt + adB(t) + dY (t). (3.44) 
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Remark 3.2. The process Q is referred to as an ROU process. The steady state distribution of Q is a 
truncated (at zero) normal variable, 


Q(o o) = Normal 



(3.45) 


whose mean is 


E[Q{ oo)] 


o V29 \ CTy/0]2J 


(3.46) 


where the hazard function h(-) is defined as the ratio 
distribution: 


h (x) 


1 — 4>(a;) ’ 


of the density and the tail of the standard normal 
for all x € TZ. (3.47) 


3.2 Preliminaries 

We conclude this section with several results that are well known in the heavy traffic literature. As such, we 
do not provide proofs. The lemmas are all similar in substance to those in Lemma 3.1 of [? ]. 

Lemma 3.3. Bounded total arrivals. For any t > 0, 


lim P (A n (t) > 2 put) = 0. 

n—f oo 

Lemma 3.4. Bounded maximum service time. For any e, K , t > 0, 

lim P ( sup vf > e/\fn\ = lim P ( sup > e/y/n \ = 0. 

n ^°° \i<Knt J n ^°° \i<Knt ) 


Lemma 3.5. Functional law of large numbers for initial service times For any e, b > 0, 


lim sup P sup 

n—foo \jik<b^/ri 


E Vi 
fl n 

2=7 + 1 r 


(k~j) 

pn 



= 0. 


Lemma 3.3 places an upper bound on the arrival process. This bound allows us to replace the number of 
arrivals in an interval with a deterministic upper bound. Likewise, Lemma 3.4 gives a uniform upper bound 
on service times that are of order n in quantity. Used in conjunction with Lemma 3.3, Lemma 3.4 places 
an upper bound on all service times during any finite interval of time. Lemma 3.5 places a bound on the 
amount by which the service times of initial jobs can differ from their expected value. 

The following four results pertain to the arrival of jobs and the arrival of potential work. The first result 
states that jobs arrive in a linear fashion. 

Lemma 3.6. Uniformly bounded fluid arrivals. For any e and t > 0, 


lim P ( sup |A"(s) - ps\ > e ) = 0. 

S<t 


The second result, based on the heavy traffic condition, is a functional law of large numbers and states 
that asymptotically, potential work arrives at rate 1, uniformly over compact intervals. 

Lemma 3.7. Uniformly bounded fluid potential workload. For any e and t > 0, 


lim P sup 


s<i 


5>r- 


> e =0. 
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The above lemma is a key component for demonstrating that server idleness is asymptotically negligible; 
see Proposition 5.9. The lemma is also used for the proof of tightness of our sequence of scaled queueing 
processes; see Proposition 5.15. However, we also need another version of the above lemma, but for short 
time intervals: 


Lemma 3.8. Net potential workload tightness. For any e and t > 0, there exists a S > 0 such that 

/ A"(«) 

lim P I sup u" — (v — u) 

n— *'°° \ u<v<t,v—u<8 . , , , 

\ ~ *=A"(u) + l 

For intuition as to why this lemma holds, consider the required server effort that all arrivals would 
contribute to the workload process if none of the customers abandoned or balked. Then center this process 
at each time t by t itself, the potential amount of time that the server could have been working had the 
server never idled. If this centered process is then scaled up by y/n and the limit is taken, the result is a 
Brownian motion. It follows that the sequence of processes is tight and this fact is used in the proof of 
Proposition 5.15 . 



4 Approximations and interpretations 

Now we use the limits of the previous section to approximate several performance metrics. Throughout this 
section, assume that we have a queueing system with arrival rate A, service rate p, balking distribution Fj,, 
abandonment distribution Fd , standard deviation of interarrival times <r Q , and standard deviation of service 
times a s . To draw connections between the formal limiting procedure with the original queueing system, we 
make the following notational conventions: 


p n = n, ,3=^—^, \ n = n + /3y/n 
y/h 


and 


a = yj(a a ■ A) 2 + (a 3 ■ pf. 

Notice that defining any two of A, p and (3 uniquely defines the third. 


4.1 Distribution of the ticket queue in steady state 

Noting the scaling Q n = \JnQ n ~ yfpQ 1 we approximate the steady state distribution of our queueing 
process using the steady state distribution of the ROU process: 

Q k Normal 0,oo^, (4.48) 

that is, the ticket queue distribution is approximated by a normal distribution with mean (A — p)/9, variance 
pa 2 /(20) and truncated to lie within [0, oo). Note the substitution: /3 ss (A — p)/yfp. 


4.2 The expected ticket queue length 

We can also approximate the expected queue length as 


where p = A /p. 
example, [? ? ? 


E(Q\ 


A — p 
0 


+ aj-h 


(1 - P) /2/f 

a V 0 


(4.49) 


If one is also interested in approximations for higher order cumulant moments, see for 
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4.3 The fraction of abandonment 


There are three approximations forwarded for the abandonment probability. To simplify notation we let 
^(0) = F^(0) and /(0) = i^(0). The first approximation takes the rate of abandonment from the queue and 
divides by the total arrival rate: 


g{0)E[Q] _ p- 1 g( 0) ag( 0) / (1 - p) [2p\ 

A p 9 p\/26p y a V 0 J 


(4.50) 


The second approach starts with computing the cumulative distribution function evaluated at the expected 
delay: 


0-1 


E 


G ( ^ 
\ x 


= E 










, d-^(O) / (1-p) [2p 

(/>-*)— + -55ttA|-t-V T 


e VW v , 

The last approach is the simplification of the first two under the assumption that p = 1: 

o'S'(O) 


«3 


2y/Tv6p 


(4.51) 


(4.52) 


4.4 The fraction of balking customers 

The balking probabilities are similar to the abandonment ones and, as such, have three versions. 


7i 


p 9 p^j29p 
72 « {p~ l)-x- + -y 



and when p = 1, we have that 


73 


°m 

2y/n 9p 


4.5 The expected number of unresolved abandoned tickets 

Given the number of unresolved tickets a fixed fraction of these are expected to be abandoned: 


E[X(t)) « Q(t) 


g(Q)Qft) 2 

2p 



9_i 0) 
2 


i?[Q(oo)] 2 . 


(4.53) 

(4.54) 


(4.55) 


(4.56) 


4.6 Interpretation 

So why should we believe that there is very little difference between the ticket queue and the standard queue 
in steady state? In the absence of balking, one would assume that the number of customers in the standard 
queue would be smaller than that in the ticket queue, as the former rids itself of customers who will not add 
to the server workload. One mechanism that reduces this difference is that the ticket queue, albeit longer, 
ultimately sees less work than its queue length would suggest. Hence it must be resolving its ticket queue 
faster than the standard queue is processing its customers. Moreover, the concentration of abandoned tickets 
is typically greater among the tickets close to the front of the queue, as these tickets have been in circulation 
the longest. But the closer the ticket is to the front, the sooner it gets resolved. The more abandoned 
tickets, the faster the server resolves such tickets. Hence, the ticket queue tends to drive itself back toward 
the standard queue status the farther away it deviates from it. 
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Adding in the balking customers further lessens the difference between the two queueing scenarios. If the 
ticket queue is longer than the standard queue then the former has more customers balking at the front end. 
To conclude, as the distance grows between the length of the ticket and standard queues, so do the forces 
that force the coupling of the two queueing models. This notion is formally expressed in Theorem 5.8. 


5 Proof of the main results 

The results that follow lead up to the proof of the main result at the conclusion of this section. Some proofs 
are delayed until the Appendix. 


5.1 Asymptotic Boundedness 

We argue first that the scaled queue length processes and the workload processes are asymptotically bounded. 

Lemma 5.1. Under (3.42), we have that for any t,rj > 0 there exists a I\ = K(i]) > 0 such that for each 
a£{S,T}, 

limsupP ( sup Qa(s) > K 

n->oo \s£[0,t] 

and 

limsupP ( sup W™(s) > I< 

n-¥ oo y«e[o,t] 

Proof. The queue length processes for the ticket and standard queues can both be bounded path-wise by a 
third queueing process that contains neither balking nor abandonment. It is standard that this third scaled 
queue length process converges to a reflected Brownian motion. It also follows that this third scaled queue 
length process exhibits the boundedness expressed in (5.57); e.g., see Lemma 3.4 of [? ]. And because this 
third process bounds the ticket and standard queueing processes for every time t, the result in (5.57) follows. 
The same arguments hold for the workload processes in (5.58) and this concludes the proof. □ 

Lemma 5.1 emphasizes the orders of magnitude of the queueing and workload processes. This lemma 
will be used frequently in conjunction with the balking and abandonment distributions to place bounds on 
abandonment and balking frequencies. 


< V 


< V- 


(5.57) 

(5.58) 


5.2 Abandonment and balking frequencies 

Next, we cover several properties of the accumulation of balking and abandonment events among the arriving 
jobs. The following lemmas, which besides Lemma 5.3, are provided without proof, use the fact that the 
derivatives of the balking and abandonment distributions exist at zero; see (2.1). The first lemma is used 
throughout this section and follows from a straightforward application of Taylor’s Expansion. 

Lemma 5.2. For any K > 0, 

Fb{K/ y/n) F d (K/ y/n) ^ 

K/y/n K/y/n 

for sufficiently large n. 

The second lemma is similar. 

Lemma 5.3. For any 6, K > 0, 

(F b ((s + 6)/y/n) - F b {s/y/n) | F d {{s + 5)/y/n) - F d {s/y/n) \ ^ OQ 

4SA ws ws ) 

for sufficiently large n. 
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The next result shows that one can choose a sufficiently small S such that, uniformly over all subintervals 
of [ 0 ,i] of size 6, the total number of jobs that arrive in any subinterval that either abandon or balk is 
arbitrarily small. The proof of this and subsequent results can be found in the Appendix. 

Proposition 5.4. For any e, rj,t > 0 and K > 0, there exists a S > 0 such that 

/ A"(s+5) \ 

limsupP [ sup E (l (bi < K/yfn) + 1 (di < K/yfn )) > eyfn I < 77 . (5.59) 

\'<t i=An(a)+1 j 

Likewise, for a sufficiently small time interval, the amount of potential workload contribution associated 
with balking or abandoning jobs arriving during the interval is smaller than order 1 /yfn. This result is a key 
element in the proof of tightness of our scaled queue length and workload processes; see Proposition 5.15. 

Proposition 5.5. For any r],t > 0 and K > 0, there exists a 6 such that 

/ \ 

limsupP ( sup "S'] v " • (1(6^ < K/yfn ) + \{d l < K/yfn )) > -^= I < 77 . (5.60) 

n ^°° \ s ^ i=Ar-(s)+ 1 v n ) 


So far our propositions have been able to replace the queueing and workload processes with upper 
bounds early in the proofs. For the following result, where we show that the centered and scaled balking 
and approximate abandonment processes converge to zero, such substitutions cannot be made immediately. 

Proposition 5.6. Centered balking and reneging processes are negligible. Under the assumptions of Theorem 
3.1, for each a £ {S, T}, and any e,i),t > 0, 


and 


lim sup P 

sup 


n—> 00 

\s 6 [ 0 ,t] 

limsupP 

sup 


n—too ' 

1 sG [ 0 ,t] 





< V 


< V- 


(5.61) 


(5.62) 


The implications here are that the balking and reneging 
respective distribution functions. 


random variables can be replaced with their 


5.3 Coupled Processes 

An interpretation of Lemma 5.1 is that the queue length is order yfn. In a model with no balking or 
abandonment, this fact would be sufficient to draw a linear relationship between the queue length and the 
workload of the form Q/p « W. In the presence of balking or abandonment, this relationship is justified in 
[? ]. The key is that the number of jobs in queue who do not contribute to the workload is negligible with 
respect to yfn. We have the same result here. 


Proposition 5.7. State space collapse. Under the conditions of (3.42), we have that for anyt,e,r] > 0 and 
each a £ { S, T}, 


lim sup P 

n—t 00 


sup 

s£[0,£] 


Q n a {s) ^ pW-(s) 


> £ 


< V- 


Next, we establish that the sequences {Qlf,n > 1} and {Qg,n > 1} converge to the same limit, if 
anything at all, as do {W^,n > 1} and { W‘f'. n > 1}. 

Theorem 5.8. Asymptotic coupling. Under the assumptions of Theorem 3.1, for any e,rj,t > 0, 


lim sup P 

n—>■ 00 


sup 

\*e[o,t] 


Qsi s ) — Qt ( s ) 



< V 


(5.63) 
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and 


lim sup P 

n—foo 


sup 

\se[o,t] 


WS(s) 


Wf(s) 


> £ \ < Tj. 


(5.64) 


Proof. We first show that (5.64) holds. Then (5.63) follows from applying the Triangle Inequality twice and 
Propositions 5.7 for both a = S and a = T. 

Fix e,rj,t > 0. Removing the absolute value signs yields 


P 


sup 

\se[o,i] 


W$(s) 


W¥(s) 



< 



sup WS(s) 
e[o ,t] 


W?(s) >e +1 


sup 

se[o,t] 


Wp(s) 


w%(s) 



We will show that 


and then 


lim sup P sup Wt(s) — Wg(s) > e I < rj/2 
n->oo \se[o,t] J 


lim sup P [ sup Wg(s)- W$(s) >e) < 77/2 

n-¥ 00 y sg [o,t] J 


(5.65) 


(5.66) 


(5.67) 


will follow by symmetry. 

There are two steps to demonstrating that (5.66) holds. First we restrict the amount by which the gap 
between Wt and Ws can grow at any instant. Then we show that once the gap is of a certain size - one 
that is smaller than e - the gap necessarily must shrink. To this end, introduce the notation for any Sq < t: 


r n (s 0 ) = inf{s G [s 0 ,t] : W^{s) - Wg(s) > e} 


and 

7 n (s 0 ) = inf{s G [s 0 ,t\ : W£(s) - W£(s) < e/2}. 

where either stopping time is equal to t if the corresponding infimum is taken over an empty set. Suppose 
Wf is greater than W§. The depletion of the former is always as fast as the that of the latter since the 
servers work at the same rate. Therefore, the gap between Wf and Wg can only increase due to jumps in 
Wf. If jumps in Wf are all strictly less than e/4, then before the gap can exceed e, it must first assume some 
value in (3e/4, e). Then, once in this interval, the process must hit e before falling below e/2. Otherwise, 
the gap must again assume a value in (3e/4, e) before it reaches e. Hence, we have that 


sup Wp(s) - Wg(s) > e < P sup W¥(s) - W^{s-) > e/4 


se[o,t] 


se[o,t] 


(5.68) 


+ P ( 3s 0 < t s.t. Wt(s 0 ) - Wg(so) G ( —,e ) and r n (s 0 ) < 7 n (so) 


3e 


For the first term on the right hand side, a jump in the ticket queue workload must be due to a large service 
time associated with an arriving job. By Lemmas 3.3 and 3.4, 


P ( sup Wt(s) — Wt(s-) > e/4 j < P (A n (t) > 2 fint) + P ( sup v n (i) > e/4 J < (5.69) 

\sG[0,t] J \i<2fint J 2 

As for the second term on the right hand side, WJf is greater than Wg throughout the interval [s 0 , min(r n (s 0 ), 7 n (so)]- 
It follows then that any new job that arrives to both queues during that interval and eventually abandons 
the standard queue must also abandon the ticket queue. So for the gap between the processes to increase 
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during this interval, it must be due to jobs that balk at the standard queue but not at the ticket queue; this 
is possible only if the ticket queue is smaller than the standard queue during this interval. It follows that 


3s 


< 


< 


0 < t s.t. W£(so) - Wg(s 0 ) £ and r ”( s o) < 7 n (so)^ 

P (3s < t s.t. W?(s) - W£(s) > | and Q£(s) - Qg(s) < o) 


«« - Wf(s) 

> 7 ) + P ( sup 

^ (s) WS(s) 

T 

4 J \ s 6[°.d 

T 



where the second inequality is a consequence of the Triangle Inequality. Applying Proposition 5.7 twice 
yields 

P ^3s 0 < t s.t. Wp(s 0 ) - Wg(s 0 ) e and T "( s o) < 7 n (so)j < (5.70) 

Hence, (5.66) follows from (5.68)—(5.70); (5.64) follows from (5.65)—(5.67) ; and (5.63) follows from Proposition 
5.7, (5.64), and the Triangle Inequality. □ 

Theorem 5.8 allows us to focus all of our efforts in proving that one of these sequences - say {Qlf,n > 1} 
converges because the other sequence is brought along with it. 

The service allocation process T n converges to the identity function. 


Proposition 5.9. Convergence of the allocation process. Under the assumptions of Theorem 3.1, for each 
a £ {S', T}, and any £,rj,t > 0, 


limsupP ( sup \T™(s) — s\ > £ j < 77 . 

n->oo yse[o,t] J 

A further implication of this result is that the sequence of idle time processes converges to zero. 


5.4 Simplifying the reneging process 

Reneging can happen only if the associated job actually joins the queue. This complication makes for 
involved expressions for the reneging processes in (3.17) and (3.19). Furthermore, notice that the expressions 
include both the queue length process as well as the workload process. The process ejj allows one to replace 
the workload process with the queue length process, to ignore whether or not the jobs have balked when 
considering whether they will abandon, and to ignore the timing of when reneging is detected by the system. 
The following proposition justifies this approximation. 

For each n > 1 and a , we define the processes = {I?° ,ra (f),t > 0}, R^ n = {R^ n (t),t > 0}, and 
R 2 U n = > o}, Where 

A n (t ) 

R°u n (t) = !(* < Q2(t?-)/M n ). (5-71) 

i =1 
A n (t) 

Ra n (t) = E !(* < TO*"-)) (5.72) 

»=1 

and 

A n (t) 

R a n {t) = E ^ > QlW-Vn ■ < W^(t?-)). (5.73) 

i=l 

Working backwards, notice that R% n is the fictitious abandonment process where the abandonment takes 
place upon arrival. Next, is the fictitious abandonment process whereby abandonment happens upon 
arrival and jobs may abandon even if they have already balked. Finally, is the process by which 

abandonment happens upon arrival, is independent of the balking process, and is based on the weighted 
queue length upon arrival instead of the workload upon arrival. The diffusion scaled analogs have the form 
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i&" = > 0}, where R^’ n (t) = (1 / y/n) R^ n (t) for each k = 0,1,2 and a £ {S, T}. Ultimately, 

we would like to replace i?” with see Proposition 5.14. The following three propositions arrive at that 
conclusion progressively. 

The first result verifies that abandonment may be treated as taking place upon arrival. 

Proposition 5.10. For any e,rj and t > 0, 

limsupP) sup (R% n ( a )-RZ(s)) > e) < rj. 
n-to o y s g[o,t] ' ' J 

An immediate consequence of the next result is that effectively, no customer could have abandoned if 
their balking random variable was small enough to cause it to balk as well. That is, the number of jobs that 
are candidates for both balking and abandonment is asymptotically negligible. 

Proposition 5.11. For any e,rj and t > 0, 


lim sup P 

n—f oo 




< v- 


Now, we show that abandonment can be based on the weighted queue length as opposed to the workload. 


Proposition 5.12. For any e, 77 and t > 0, 


lim sup P 

n—too 


sup 

K n (s)-K n (s) 

y « 6 [o,t] 



> e 


< 77- 


Proposition 5.13. Reneging among initial jobs is negligible. For every rj > 0 and each a £ {S, T}, there 
exists an L > 0 such that, under the assumptions of Theorem 3.1 , 

limsupP (Ra(t) > L\ < r}. 

n—> oo ' ' 


Proposition 5.14. Reneging effectively ignores balking and takes place upon arrival. For each a £ {S,T}, 
and under the assumptions of Theorem 3.1, 

-> 0 


in probability as n —► oo. 


Proof. Notice that for any t > 0 and a £ {S', T} that 


m = 



(KW - ;%"(*)) + (K n (t) K n (t)) + (K n (t) - . 


The result is an immediate consequence of Propositions 5.10, 5.11, 5.12, and 5.13. 


□ 


5.5 Tightness 

Next, we argue that the scaled processes are tight. This result is a key step in proving Proposition 5.17, 
which in turn is used to extract the restorative drift of the limiting diffusion process. 

Proposition 5.15. Tightness of the scaled processes. The processes > 1} and { Wf , n > 1} are tight. 

Proof. By Theorem 13.2 of [? ], tightness follows from Lemma 5.1 and the fact that for any e > 0, we have 
that for either a £ {S,T}, 


lim lim sup P 

^0 n—too 


sup 

\u,v(z [0, £],?;— u<8 


Qliv)-Q n a iu) 



= 0 


(5.74) 
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and 


lim lim sup P sup 

<5->-0 n_>oo \ U,ve[0,t],v-u<6 


- WZ(u) 


> e = 0. 


(5.75) 


We will show that (5.75) holds. Then (5.74) follows from Theorem 5.7 and (5.75). This will conclude the 
proof. 

Fix a £ {S', T} and arbitrary constants e,r] > 0. We will show that there exists a do > 0 such that for 
any S < d 0 , 


sup \W2(v)-WZ(u)\>—= )< V 

i u,v£[0,t],v—u<5 V ^ / 


(5.76) 


for sufficiently large n. Changes in the workload are due to the arrival of work and the processing of work. 
That is, for any v > u > 0, 

A n (v) 

WS(v) W2(u) = E v i ■ ^ > QlW-)/n ■ 1(* > - (T a » - T™(u )). (5.77) 

i=A n (u )+1 

By (3.23), the change in the service allocation process can be written in terms of the increase in the idle 
process, which in turn can be bounded by the shortfall in the arrival of workload relative to the potential 
server effort. For any v > u > 0, 


i2(«)-i2(u) = (v-u)-(T2(v)-T2(u)) 

A n (s) 


(5.78) 


< max 0,- inf E < • 1(6, > ' 1 (* > W^-)) - (* - «) 


It follows from (5.77) and (5.78) that for any 0 < 5 < do, 


sup \W™(v) - W™(u)\ < 2 sup 


u,v€ [0, t] ,v—u<.8 


u,v£ [ 0 , £],?;— u<8 


A n {v) 

E «?-(«-«) 

i=A n (u)-\-1 


(5.79) 


A n (n+<5) 


+2 sup E v i U& ^ QWi-)/»") + Udi < TO*"-))) 
«e[o,t] i=A „( u ) +1 


< 2 sup 

IL,V (z[0, t],v — u <8 0 


A n (v) 

E «? -(v-u) 

i=A n (n) + l 


A n (n+^o) 

+2 sup E < < Q n a.W-)/n + 1(* < TO*?-))). 

1=A»(«)+1 
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Finally, 


sup \WZ{v)-WZ{u)\> — 

Lt,i>£[0,t],'U— u<5 V 7 , 


< 


2 sup 

u,v£ [0,t],u— U<Sq 


A n ( v) 

Y v? ~(v-u) 

i=A n (u)-\-1 
A n (u+<5q) 


> 


2 y/n 


< 


+P 2 sup Y v i Wk ^ QWi~)/l* n ) + *(* < WS (*"-))) > 77 

A" 7) 

^ ( v - u ) > 

i=A n (u)+1 
A n (n+<5o) 


2 sup 

ii,'PE[0,£] , 1 ;—n<«5o 


2^/n 


sup Q™(s)/p n > K/y/n 

s<t+8o 


( \ UJ \ / \ 

2 sup V vj* ( 1 ( 6 , < Jf/Vn) + l(dj < K/yfn)) > -^= + P ( sup W"(s) > K/y/n ) . 

^i=AHu) + l 2 ^J W*+«o / 

By Lemma 3.8, Lemma 5.1, Proposition 5.5, and (5.79), we can choose a, K > 0 and a <5o > 0 such that 
(5.76) follows easily. 

□ 


5.6 Convergence to Diffusion Processes 

In this subsection, we pull together the last ingredients for our diffusion process. We start with the driving 
Browning motions from our centered and scaled interarrival and service time processes. We derive the 
restorative drift term of our limiting diffusion process from the derivatives of the balking and abandonment 
distributions evaluated at zero. Finally, we complete the proof of our main result. 

Let B a and B s be independent, standard Brownian motions; i.e., B a (0) = B s ( 0) = 0 and the processes 
have zero drift and unitary infinitesimal variance. The diffusion scaled arrival processes and service comple¬ 
tion processes converge to scaled versions of these Brownian motions. These are standard results - see, for 
example, [? ]. 

Proposition 5.16. Under the assumptions of Theorem 3.1, 

A n => pa a B a and => pa s B s 


as n —> oo. 

The process d” swaps the sum of the abandonment and balking distributions, evaluated at the scaled 
queue length at the times of arrivals, with a smooth function involving the derivatives of the distribution 
functions evaluated at zero and then multiplied by the scaled queue lengths. The following proposition 
justifies this step and is proven in the Appendix. 

Proposition 5.17. For each a € {S, T}, and under the assumptions of Theorem 3.1, 

^ 0 


in probability as n —► oo. 

We conclude with a proof of the main result. 

Proof OF Theorem 3.1. Fix a £ {S', T}. Recall the expression for the scaled queue length and 
idleness process given in (3.39). The elements of the process X™ (see (3.40)) converge to either zero, a drift, 
or a Brownian motion. By Lemma 3.6 and Proposition 5.6, we have that 

M{ 0 od n AO and 
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as n —> oo. Furthermore, Propositions 5.14 and 5.17 inform us that 


eJJ —0 and S™ —> 0, 


respectively, as n —> oo. Proposition 5.16 provides the convergence of the scaled and centered arrival and 
potential departure processes. Coupled with the service allocation process, we have 

A n + S% o TZ ->• p<j a B a + na b B b = oB. 


Finally, (3.36) provides the drift term 


(A 11 - n n ) 




as n —> oo. Hence, 

X%^X. 


From (3.42) we have that 


(Q n (0),XZ)^(Q 0 ,X), 


and hence, by the Continuous Mapping Theorem, 


as n —>• oo, 


(<&,*?) = (4> e ,4/ e )(Q"(0),A-) =► ($ 9 ,* 9 )(Q(0),X) = (Q,Y), as noo, 
where recall Q(0) is equal in distribution to Qo- 

As a was chosen arbitrarily, this limit holds for both the ticket queue- and standard queue-related 
processes; both scaled queue length processes converge to the same limit in distribution. By (5.63) of 
Theorem 5.8, these processes are coupled: 


((<& Yg), ( Qt , Yt)) =► ((Q, Y), (Q, Y)), 


as n —> oo, 


Finally, Proposition 5.7 demonstrates that asymptotically the scaled workload and scaled queue lengths are 
scalar multiples of each other and thus converge together to scalar multiples of the same process. Hence 
(3.43) holds. This concludes the proof. □ 


6 Numerical results 

We have tested the heuristics forwarded in this paper extensively. Tables capturing these tests are given in 
the Appendix. We describe some of the highlights and trends in the results here. In Tables 1 - 4, we see 
that our approximations generated from the heavy traffic limit theorems are quite good, especially when the 
arrival and service rates are large i.e. /i ss 100. In Table 1, we simulate a scenario where all of the parameters 
are generated from exponential distributions. We see that the approximations are good and accurate, except 
when the rates are small and the drift parameter f} is very negative. 

In Table 2, we simulate a non-Markovian model where the service rate follows a lognormal distirbution 
and the balking and reneging distributions follow a uniform distribution. Unlike the exponential distribution, 
the uniform distribution is bounded and we see that the simulated values in Table 1 are very similar to the 
values of Table 2 even though they differ in the types of distributions that generate the queueing dynamics. 

In Table 3 we simulate a scenario where all of the parameters are generated from exponential distributions. 
This table is different than Table 1 since we use different balking and reneging parameter values. Once 
again we see that the approximations are good and accurate, except when the rates are small and the drift 
parameter /3 is very negative. In Table 4 we simulate a non-Markovian model where the service rate follows a 
lognormal distirbution and the balking and reneging distributions follow a uniform distribution. Once again 
the simulated results of Table 3 are very similar to Table 4 even though the distributions are not Markovian 
in Table 4. 

In addition to our heavy traffic approximations accurately estimating the performance measures, we also 
notice that in every simulation the ticket queue is larger than the standard queue. This fact is irrespective 
of the distributions that are used to generate the queueing dynamics. Moreover, the fraction of balking in 
the ticket queue is also larger than the standard queue. Even though these two processes converge as the 


21 



rates tend to infinity, when the rates are finite, the ticket queue is perceived as being larger, which causes 
more people to balk from the system. This difference in queue length is also a function of the parameter fi. 
When the ft parameter is larger and positive the difference is larger than when /3 is negative. Thus, when 
the rates are not infinite, the ticket queue exhibits interesting behavior that is to be expected. 

7 Conclusions and extensions 

In this paper we studied the dynamics of a critically loaded queueing system subject to customers who might 
either balk because the line is too long or abandon the queue from waiting too long in-queue. We consider 
two types of abandonment protocols. In the conventional approach to capturing abandonment, customers 
immediately leave the system when their patience time has been exceeded by their time waiting for service; 
we refer to the model in this setting as the standard queue. In the standard queue, everyone is aware of 
an abandoning customer’s departure at the time of abandonment. We compare the standard queue to the 
ticket queue. In the ticket queue, customers whose patience has run out leave the queue in an unnoticed 
fashion. Their departure is only detected when their hypothetical service time would have begun. The paper 
is complementary to the study of [? ] who study a heavily loaded system with impatient customers; in 
comparison, our customers are relatively patient. Our method of analysis is also fundamentally different 
from that of [? ]. 

We prove a heavy traffic limit theorem for the diffusion scaled queue length and workload processes. A 
key result in the theorem is that the standard and ticket queues are asymptotically coupled under diffusion 
scaling. The managerial interpretation is that regardless of how you implement your queue - whether with a 
physical line, which is best modeled with a standard queue, or as a ticket queue - the dynamics of the systems 
will not differ by much. In addition to this insight based on the sensitivity of the diffusion scaling, we provide 
some heuristics for calculating certain performance metrics of operational importance. These heuristics are 
beyond the sensitivity of diffusion scaling. Nonetheless, we assess the accuracy of these heuristics through 
simulation. We find that in a broad range of parameter and distributional settings, the heuristics perform 
well. 


8 Appendix 

The Appendix is split into two parts. The first half provides proofs of results stated earlier in the paper. 
The second half provides an extensive collection of numerical examples. 

8.1 Proofs 

This Appendix contains many of the proofs from Section 5. We start with a proof of Proposition 5.4. 

Proof of Proposition 5.4. Fix e,r],t,K > 0 and set 6 = A ^ Ke . By Lemma 5.2 and our choice of 8 , 
we have that _ 

/m8 (Fb(Ky/n) + Fd(K/y/n)) < 2 ^\fnK98 < 
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for sufficiently large n. It follows that, by Lemmas 3.3 and 3.6, 


A 71 0+5) 

sup ^2 (l (h < K/y/n) + 1 (di < K/y/n )) > ey/n 


i=A"(s)+l 


< 


j+[2p,n6] 

sup V' (l (bi < K/y/n) + 1 (di < K/y/n)) > ey/n 

0<j< \2y,nt\ i= j +1 


+P ( A n (t) > 2 /mt) + P ( sup {A n (s + S) — A n (s)) > 2 /m5 

S<t 


< / + P 2 sup 


0<t<L2t/5J + l i=jL2M 


+P ( 2 sup 

0<j< \2t/5\ +1 


+P ( 2 sup 

0<j< \2t/S\ +1 


0+l)L2/ira5J 

^ (1 (bi < K/y/n) - F/K/yfn )) > ^ 

2fin8\-\-l J 

Cj+i)L 2 a*toJJ /—\ 

^ (1 (di < K/y/n) - F/K/y/n)) > ^ 

i—j [2+n5J +1 / 

(j+l)l2/j,nSi ,—\ 

^ F/K/y/n) + F d (K/y/n) > ^ 

i=j\2^nS\+l ) 


for sufficiently large n. 

We will show that, for sufficiently large n, 


0'+l)l+«<5J 

2 sup y2 (l(&* ^ K/y/n) - F b (K/y/n )) 

0<j<L2i/5J + l i=j ^ nS ^ + 1 


> 


ey/n l n 
6 < 6 ’ 


2 sup 


o<i<L2t/*j+i. =J>n5J+1 
and by symmetry it will also follow that 


^ F b (K/ y/n) + F/K/yfn) 


> 


Ey/n l n 

6 < 6 


/ C?'+1)Lm™‘5J ,—\ 

2 sup Y' (l (di < K/yfn) - F/K/yfri]) > C ” 

V °- j -t 2t/Si+1 i=jlfi.nS\+l 6 J 


Ey/n l ?) 

< +• 


Consider (8.81). By Kolmogorov’s inequality (see, e.g., [? ]), 

(.7+1) La 4 ™ 5 ! 


2 sup 
0<J<|2i/«J+l j=J - L/i . 


1)L/Jn5j \ 

Y2 (l(6i < K/y/n) - F/K/y/n)) > ^ 
L+niJ+l J 


< 


2 1 


144 


E 


L+n5J 

Y2 ^ K /Vn) - F/K/y/K)) 

2=1 


(8.80) 


□ 


(8.81) 


(8.82) 


(8.83) 
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The indicators above are independent, so the cross terms all have expectations of zero. Hence 


sup 


0<j<L2/mt/SJ+l i=jL#m(5J+1 


j— 

E (U 6 * ^ A '/v^) - F b (K/Vn)) > 


< 


< 


E 


\_finS\ 


£ (l(^< K/y/n)-F b {K/yfn)Y 


< 'i 
- 6 

for sufficiently large n. So (8.81) holds, as does (8.83) by symmetry. The result (5.59) follows from (8.80)- 
(8.83). 

Next we prove Proposition 5.5. 

Proof OF Proposition 5.5. We will alter (5.60) to an expression that, without loss of generality, 
excludes the abandonment random variables. We will show that for any 77 , t > 0 and K > 0, there exists an 
S > 0 such that 

/ A"(s+<5) \ 

P SU P E V i ’ 1 ( & * ^ K N n ) > 7 = < ( 8 - 84 ) 

V s -* i=A"(s)+l J 

for sufficiently large n. Fix rj,t and K and set S = 12 e K g and notice that by Lemma 5.2, 


i^ n ) F b {K/y/n) < 4S ^ < — 

VA 4 / V 71 3 V n 

for sufficiently large n. Proceeding in a manner similar to that of the proof of Proposition 5.4, it follows by 
Lemmas 3.3 and 3.6 that 


a"(h-<5) 

sup E v i ' l (Pi < K/y/n) > ~^= 
• s - 4 »=A"(s)+l * 12 


< 


j+limSj 

sup E v ? ' (!( & i < K /Vn) + 1 (di < K/y/n )) > ~^= 

0<j< [2fint\ i—j + i V n 


+P ( A n (t ) > 2 fj,nt) + P ^sup (^"( s + 5) — -A n (s)) > 2 finSj 


< - 


• j+\_ixn8\ , . 

+ P sup E ( v i ~ — ) ■ Y b i ^ K /Vn) > E- 
\o<j<\2^nt\ .frf. V A 4 "/ Wn 


y0<j<[2g,nt} i—j + i 

j+l/mSj 


+p ( sup E 4- (Wi - F/'M - F b (I</y/n )) > —E 

\o<?<L2/mtj i ±~j_ 1 A 4 V 74 


(8.85) 
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for sufficiently large n. Consider the third term on the far right hand side. By Kolmogorov’s inequality [? ], 


( j+[fmS\ 
sup 

0<j<L2/intJ i—j + i 


■ (10 i < K/y/n) - F b (K/y/n)) > 



< 

< 

< 

< 

< 



sup 

0<7<L2V<5J+1 




i n ^) F b( K /Vn) 

8K6\ 


( 8 . 86 ) 


for sufficiently large n. As for the second term on the last right hand side of (8.85), applying Kolmogorov’s 
inequality a second time yields 



where the last inequality follows because of the independence of the cross terms makes their expectations 
zero. Further the service times are independent of the balking random variables. Hence, 



< 


< 



(2 1 + 25) 



(nfiS) F b (K/y/n) 

f2K6\ 

M / V Vn ) 


< 


V 

3 


( 8 . 88 ) 


for sufficiently large n. The result (8.84) follows from (8.85)-(8.88). 

□ 

Proof of Proposition 5.6. It suffices to prove that the centered and scaled abandonment process 
is asymptotically negligible. The analogous property for the balking process can be proved in an identical 
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fashion. Fix e,r],t > 0. By Kolmogorov’s inequality we have that 


sup M/ a (s) > e 

sS[0,t] ) 


SUP M£ q (s) > yfc. 

se[o,t] 


< 


ru? E 


(■ 


lnt\ 


E m < Q2(*?-)//o - F d {Qiw-)/n n )} 


Now by Burkholder’s inequality and bounds for indicator functions, there exists a c > 0 such that 


P sup MS ia {a) > e 
\*e[o,t] > 


< 



ne 2 


c „ 

< 

— 


ne 2 


2 

< 

—P 


e 2 

< 

V 


LntJ 


E t 1 ^ < - F d {QlW-)/n n )f 


lnt\ 


E m < QlW-)/n n ) + F d [QlW-)/n n )} 


di < ma xQ™(s)/n n 

S<t 


for sufficiently large n. The last inequality follows from Lemma 5.1. This completes the proof. □ 

Proof of Proposition 5.7. Fix a £ {S,T} and t > 0. For each s > 0, let Q2( s ) denote the difference 
between the index of the last arriving job and the index of the job currently in service, plus any of the initial 
jobs present at time zero that remain in the system. The jobs present at time zero do not have indices. 
Note that, for each s > 0, Qa( s ) > Qa( s )> f° r both ticket and standard queues. The process Q ” ignores 
the balking and abandonment that has taken place since the arrival of the job currently in service. One 
can think of the process as progressing in a manner similar to a ticket queue for which, in addition to the 
abandoned tickets, balking is not accounted for until service would have begun for the departing job. We 
can bound the difference between the processes. Fix an arbitrary 77, e > 0 and choose an L > 0 such that by 
Lemma 5.1, 


P( sup (Q n a (s) - Q n a (s)) 
\se[o,t] v y 


> 


< P ( sup Q”(s) > Lsfn\ + P [ sup !F”(.s) > L/y/n ) 

\se[o,t] ) \«e[o,t] / 

/ j + l(L+e/3)Vn\ \ 

+ P sup E (1 (bi < Ly/n/ 1 / 1 ) + 1 (di < L/y/n)) > E™ 

i=J+1 3 J 

f 3+ L(L+e/3)V«J 

< E P sup E (l (bt < L/y/n ) + 1 (di < L/y/n)) > 

5 i=j+1 3 


(8.89) 


for sufficiently large n. To appreciate the above inequality, note that Q 7 / is an inflated version of Q The 
jobs in the former not accounted for in the latter must have abandoned or balked, or will have eventually 
abandoned. So if there exists an s £ [0, t] such that Qa( s ) ~ Qa( s ) exceeds (L + e/3 )y/n and Q™(u) < Ly/n 
for all u £ [0,t], then there must be at least ty/n /3 abandoned or balked jobs within some (L + e/3 )y/n 
consecutively arriving jobs. We place upper bounds on the queue length and workload and this makes our 
abandonment and balking indicators i.i.d. random variables. 
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For any S > 0 it is true that (L + e/3 )y/n < fiSn for sufficiently large n. Hence, by Proposition 5.4, 


j+l(L+e/3)^n\ r- 

sup Y] (l (bi < L/y/n) + 1 (di < L/y/n)) > 

l=j+ l 3 


< 


< 


( j+[n5nj 

SU P Y (! ( bi - L /^n) + 1 ( d i ^ L /Vn)) > ^ 7 - 

i=j + 1 3 j 


A n (s+28) 


-CT. .— 

sup Y (l (bi < L/y/n) + 1 (di < L/y/n)) > 


s -* »=A"(s) +1 


+P ( inf (A"(s + 25) - A n (s)) < fiSn 


< i 


(8.90) 


Additionally, define for each s > 0 the quantity W/(s) that tracks the service times of the jobs associated 
with Qa(s). This process has the entire service time of the job currently in service and the service times of all 
jobs that arrive after the arrival time of the job in service; jobs that abandon or balk contribute the process 
W£ nonetheless. Just as with and its augmented version Q™, it is also the case that W/(s) > W£(s) 
for each s > 0. Unlike the process lU”, the augmented process W™ experiences both upward and downward 
jumps. Upward jumps are the size of the would-be service time of each arriving job, even those that balk or 
abandon, and occur at the time of arrival of the corresponding job. The downward jumps occur at service 
completion times. The downward jump size is equal to the service time of the job that was in service plus 
the would-be service times of jobs that arrived between the arrival time of the job just served and the arrival 
time of the next job to be served. Just as we did in (8.89) for the queue lengths, we can bound the difference 
between these proceses. By Lemmas 3.3, 3.4, and 5.1, and Equation (8.89), 


P ( fi n sup fe(s) - 


> 


< P U n • sup < > 
\ i<A n (t ) 


£y/n 


sup Q”(s) > [L + e/3)y/n 
^se[o,£] ) 


(8.91) 


+P sup Q™(s)/fj, n > L/y/n ] +P sup W/(s) > L/y/n 


se[o,t] 


se[o,t] 


j+ LU+e/3) 

+ P | n' 1 ■ sup V' v/ ■ (l ( bi < L/y/n) + 1 (di < L/y/n)) > 

j<A n (t) i=j+1 

j+ [(L+e/3)y/nj 


< - + P | M ' SUp 


Y v i • (! ( b i < L /Vn) + 1 (di< L/y/n)) > 
o<AHt) l=J+1 b 


for sufficiently large n. As with (8.89), we replace (L + e)y/n with a bigger quantity nSn (provided n is 
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sufficiently large), where, by Proposition 5.5, <5 > 0 is chosen such that 


i+L(L+e/3)^/nJ 


A i n sup v™ ■ (l (bi < L/y/n) + 1 (di < L/y/n)) > 

j<A n (t) i=i+1 


i+L/i<5nJ 


< P U" sup y' v" ■ (1 (bi < L/y/n) + 1 (di< L/y/n)) > 


T"(s+ 25) 


< P|/i n sup y2 v™ ■ (l (bi < L/y/n) + 1 (di < L/y/n)) > 


i=A"(s) + l 


+P inf ( A n (s + 25) - A n (s)) < fi5n 

' S<t 


ty/n 

6 

e^/n 

nr 

£y/n 

nr 


< -7 


for sufficiently large n. 

Lastly, note that by the functional weak law of large numbers, 


P sup 
V sG[0,t] 


Ql(s)-n n WZ(s) 


> 


S\j n \ n 

nr) < r 


for sufficiently large n. 

It now follows from the Triangle Inequality and (8.89)-(8.93) that 


P 


sup 

\s6[0,t] 


Ql(s) - nWZ(s) 



P ( sup 

Q n a(s)-Ql(s) 

£y/n \ 
> \ 

+ P [ fj, n sup 

\ sG[ 0 , 1 ] 


3 ) 

\ sG[ 0 ,t] 


WZ(s) - W2(s) 


> 



+P 


sup 

*e[o ,t] 


Q n a (s)-n n WZ(s) 


> 


Ey/n 


< V 


for sufficiently large n. This completes the proof. 

Proof of Proposition 5.9. Fix E,r],t > 0 and a. The server cannot work faster than rate 
follows that for each s <t, 

Ta(s) < s. 

Using (3.22) we can provide a lower bound on the service allocation process for any s > 0, 

t-( S ) > -wz( S ) + £ v? - yy v? (i ^ < Q n a (ti-)/n n ) +1 (^ < wz(ti-))). 

i—1 i— 1 

It follows from 8.94, and (8.95) that 


P 


sup 

\se[o,i] 


i to 



< 



sup T” (s) < S — £ 
e[o ,t] 



sup W£(s) > 
e[o A 


/ 

A n (s) 

sup 

£ <-* 

\ «G[0,t] 

i=l 


> 


E 

3 


(A n {t) 

+p ( yy < ■ m < Qi(u-)fn n ) + m < wto))) > e - 


(8.92) 


(8.93) 


O. 

one. It 

(8.94) 

(8.95) 

(8.96) 


28 





















From Lemma 3.7 and Lemma 5.1 we can bound the first two terms on the right hand side, 


sup W™{s) > - + P sup 


se[o,t] 


se[o,t] 


A"( S ) 

E»?-» 

4=1 


> „ I < 


(8.97) 


For the third term, by Lemma 3.3 and Proposition 5.5, 

( A "(t) 

E «" ■ ^ < Qa(ti-)/n n ) + l(di < wsiu-))) > f 


< 


< 




2 lint 


’ (A n (t) > 2 tmt) + P E ' (!& < Q2(U-)/n n ) + m < W2(ti-))) > 3 

V 2=1 

/2 [int , \ \ 

' ^ A (i(&i < Q n a (U-)/» n ) + 1(di < w^(ti-))) > " 


(8.98) 

(8.99) 


E 

\ 4=1 

/2fint \ 

( E — • (!(^ < «H/m") - F b {Ql{U-)/n + 1 (di < ^(^-)) - ^(W^-))) > 

/ 2 lint \ 

( E — ' (^(QS(*i-)/M") + ^(W^(*i-))) > |J +P(A"(t) > 2 M rrf) 


< + ^ 
“ 8 8 8 8 


< d 


The result follows from (8.95) - (8.98) and a modification of the proofs of Propositions 5.4 - 5.6. □ 

Proof of Proposition 5.10. Fix e,77 and t > 0. First notice that for any n > 0, s > 0, and 
a € {S.T}, we have that R% n (s) — i?”(s) > 0. It is instructive to expand this difference for each of the a 
values. For the standard queue, 

A n (s) 

R^ n (s) - R n s (s) = E > QsW-)/n n ) ■ l(di < W£(t?)) • 1(4 >»-*?) 

2=1 

A n (s) 

< E ^ < W sW)) • W(*?) > * - *?)• 

4=1 

Note that we have eliminated the indicator associated with balking. Moreover, for an abandoning customer 
the workload upon arrival must exceed the patience quantity. Hence we can replace the patience quantity 
in the last of the indicators with the workload upon arrival. Now we consider the ticket queue: 

A n (s) 

- RZ(s) = E W* > Q n At?~)/v n ) ■ Kdi < W?(i?)) • W(C) > s - i?) 

2=1 

A n (s) 

< E ^ E W r(t ")) • > a - *?)■ 

4=1 

Both the standard and the ticket queue have the same bounds: 


P sup 

VsG[0,t] 


R% n (s)-R n a (s) 


A n (s) 


>e < 


sup E U* < TO*"))' 1 (W£(t?) > s - i?) > \/ne 

«e[o,t] i= i 
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for each a £ {S, T}. For the remainder of the proof, fix a. Next we replace the workload quantities with an 
upper bound: 


P sup 
\se[o,t] 


R 2 a n (s) - R] 


:(s) >6^ < p( sup y' l(di < K/y/n) ■ \(K/y/n > s — t”) > y/ne | 

j \-6[o.t] iri ) 


+P sup W^(s) > K/y/n . 

\se[o,t] J 

Notice that in the first term on the righthand side above, the only jobs that contribute positively to the 
summation are those jobs i whose arrival time is after s — K/y/n; that t™ > s — K/y/n. By Lemma 5.1, 


P sup 
\se[o,t] 


R^ n (s)-RZ(s) 


> e ) < /; + P | sup 
se[o,t] 


A n {s) 

Y. 1 {di < K/y/n) > 

i=A rl ((s-K/y/n)+) 


for sufficiently large n. For any arbitrarily chosen <5 > 0 it is true that 5 > K/y/n for sufficiently large n. It 
follows then by Proposition 5.4 that we can choose a 6 > 0 so that 


P sup 
V s 6 [ 0 ,t] 


R 2 a n (s)~K(s) 


> e I < ^ + P I sup 


A n (s) 

'y 1 (di < K/y/n ) > y/ne 


[°,t] i=A n.(( s _ K /yfii) + ) 

( 

< ^ + P I sup E 1 (di < K/y/n) > y/ne 

\j>e[0,t] i=A n(( s _ 5 )+) 

< V 


for sufficiently large n. This completes the proof. □ 

Proof of Proposition 5.11. Fix e, 77 , t > 0 and a £ {S', T}. Notice that R^; n — R/; n is nondecreasing. 
Hence, replacing the workload and queue length with an upper bound yields 


P sup 

\s 6 [ 0 ,i] 




R* n (s) - i% n (s) > H < P E W* ^ Q2(ti-)/» n ) ■ *K < W?(t?-)) > Vne 


[2fintj 


— P ( E IK — K /Vn) • 1 (di < K/y/n) > /ne I +P(A n (t) > 2 /mt) 


+P sup Ql(s)/n n > K/yfc +P sup W2(s) > K/yfr . 


,«e[o,t] 1 \se[o,q 

Choose a K > 0 such that by Lemma 3.3 and (5.57) and (5.58) of Lemma 5.1, 

/ \ /\2^nt\ \ 

P ( su P] \Ryis) - 7% n (s)| > A < \ +P ( E !( 6 < ^ A 7v^) • IK < K/y/n) > /Tie J (8.100) 

As for the second term in the right hand side above, we resort to adding and subtracting the mean of 
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each summand, 


I ^2 l(^i < K/y/n) ■ l(di < K/y/n) > y/ne 


\ *=1 


/ [ 2 /intJ /— \ 

< P I ^ E [1(6, < K/y/n) ■ 1 {di < K/y/n)\ > ^ I 

+P ( ^ (1(6, < K/y/n) • 1 {di < K/y/n) - E [1(6* < K/y/n) ■ 1 (d, < K/y/n)]) > ^ I 


< P ( (2 fint)F h {K / ^/n)F d {K / ^/n) > 


/ne 


( [2fint] j— 

Y, ( X K ^ A '/Vn) • IK < K/y/n) - F b {I</y/n)F d (K/y/n)) > ^ 


By Lemma 5.2 we can bound the first term on the right hand side: 


(2 nnt)F b {K/Vn)F d {K/Vn) < 


( 8 . 101 ) 


( 8 . 102 ) 


for sufficiently large n. As for the second term, by Chebyshev’s Inequality and (8.102) 


[ 2/int J 


y" (1(6,: < K/y/n) ■ 1 (di < K/y/n) - F b {K/y/n)F d (K/y/n )) 


> 


< 


< 


rE 


ne* 


[2fint\ \ 

^ (1(6, < K/y/n) • 1 (^ < AV^n) - F b {K/y/n)F d {K/y/n)) 


-E 


2—1 






^ (1(6, < K/y/n) • 1 (^ < K/y/n) - F b (K/y/n)F d (K/y/n ))“ 


2=1 


< -F(2fj,nt)F b (K/y/n)F d (K/y/n) 

ne z 

< \ (8-103) 

The result follows from (8.100) - (8.103). □ 

Proof of Proposition 5.12. Fix e,rj,t > 0, a £ {S',! 1 }, and set (5 = e/(l6/A). By Lemma 5.1 and 
Proposition 5.7, respectively, we can choose a I\ > 0 such that for sufficiently large n, 


and 


sup W/{s) > K/y/n \ < min f^, 
se[o,t] J V 4 128/it/ 


sup \QZ(8)/n n -W2(8)\>6/y/a] <1- 
se[o,t] J ° 


(8.104) 


(8.105) 
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Fix such a I\. We start by replacing the weighted queue length with the workload process: 


P sup 
V «e[o,t] 


R 0 a n (s)-R 1 o ; n (s) 


> e 


= P ( sup 
se[o,t] 

U™(t) 


A"(s) 


E !(* < ~ m < WZW-)) 


> \Jne 


< P E W e l W ZW~) - t/Vn, WZW-) + 8/Vn}) > Vne 


+P sup \Q n a {s)/^ - WZ{s )| > 5/Jn 
\ s e[o,t] J 


(8.106) 


The second term of the far right hand side is handled by (8.105). As for the first term, we construct a 
Martingale. By Lemma 3.3, 

(A "W 

E m S [WZW-) - 5/^ WZ(t?-) + 5/V ^]) > \fne 


i=l 


/ \2^int\ 


< - 


7 {+ p [ E E[l(d i €[W?(t?-)-J/VS,TOt?-)+VVnl)|J?-i] 


> 


/ne 


\ i =i 
' [2fintj 


+P E (l(dielW2(t?-)-6/Vn,WS(f?-) + 6/Vn}) 


2—1 


-E[l(d i GTOtM-V^TOt?-) + <y/Vn|)|^r-i]) > 


/ne 


E {FdiWSW-) + 6/Vn) - ^(TO?-) - 5/y/n)) 

( [2/j.ntJ 

E e [W2(*?-) ~ W”(i?-) + <f/Vn|) 

- (^(W^(2?-) + S/Vn) - Pd(W^(t?-) - 8/y/n))) > 


> 


/ne 

2 


(8.107) 


As for the second term on the far right hand side of (8.107), notice that each summand contributes an 
amount equal to the increase of the abandonment distribution function over an interval of length 25. By 
(8.104) and Lemma 5.3, 


[2fint\ 


E (^(TO?-) + 5/Vn) - F d (W2(t?-) - 5/Vn ])) 


> 


/ne 


i—1 


< 


sup W n (s ) > K I + 1 I 2\xnt sup (F d ((s + 2 5)/y/n) — F d {s/y/n)) > 


s€[0,t] 


«e[o,ir] 


< i 


(8.108) 

for sufficiently large n. Now consider the third term on the right hand side of (8.107). The following steps 
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are similar to those in the proof of Proposition 5.6. By (8.104) and Lemma 5.2, 


E (!(* e (*?") - */V», + 5/Vn]) 

- (Fd(W2(t?-) + S/M - - 8/M)) > 


/ne 


< 


-E 


ne* 


\2fint\ 


E (l(di e FO*?-) - 8/y/n, WZW-) + S/M) 

2=1 

- (Fd(WZ(t?-) + S/M - F d {W2(t?~) - 8/M))) 2 


< —yE 

ne z 


[2/intJ 


E (!(* e [^« (*?-) - VV7 + VVn]) 


- (F d (W2(t?~) + S/M - F d (W2W~) - 5/M)Y 


< 


< 


rE 


ne* 


16/it 


[ 2fint J 

E ^(^ n (C-) + <5/v^) 


2=1 


sup W%{s) > K/y/n \ + F d ((R' + 8)/M) 

se[o,t] ) 


< 4 


V 

4’ 


(8.109) 


for sufficiently large n greater than ^ — ^ _ 

The result follows from (8.105) - (8.109). □ 

Proof of Proposition 5.13. Fix rj,t > 0 and a £ {5, T}. By Lemma 5.1 there exists a > 1 such 
that, for sufficiently large n, 


P 



sup W£(s) > 
«e[o,t] 




We use these constants to replace the workload and queue length: 


P(^(t)>i) < 


< 


< 


4r(o) 

E m < Mi) > l 

2=1 

/ 1 K V™1 


( 8 . 110 ) 


/ L \ _ / TV- \ 

I E X (7 < A 7v^) > L I + P (Q n (0) > KM + p ^(0) > 




p 


LA'\AIJ 

E 1(7 < K/M) > L 


1 

2 


for sufficiently large n. 

Recall that the residual deadlines of the initial jobs may have different distributions, f/i.ii but those 
distributions have a common bound near the origin. Namely, there exists an / > 1 and an ho £ (0,1//) 
such that Fd,i(h) < hf for each h < ho- Given such an ho and /, set L so that L > max.(2fK 2 ,4fK/^/fj). 
We replace the initial jobs’ residual deadlines with random variables which are i.i.d. In particular, on the 
same probability space, each di is replaced by di- Whenever di < ho then </></. Moreover, the common 
distribution of each di is 


F(x) 


xf, x < ho 
1 , x > ho- 
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When K/yJn < ho we have that 


(Ky/n) F{K/y/n) = fK 2 < L/2 
so that by Kolmogorov’s inequality [? ], 


/ L Ky/n\ 




I ^ 1 (di < K/y/n) > L J < P I ^ 1 (di < K/y/n) > L 


y *=i 


( 8 . 111 ) 


< 


< 


y *=1 

/L*V™I 

E (l(^; < K/y/n) 


y *=1 

4 

ZS 


n)1> § 


(/fv/n)E[(l(di < K/y/n) 

< {/) Mf(*/V5) 


< 1 

for sufficiently large n. The result follows from (8.110) and (8.111). 

Proof of Proposition 5.17. We can break 5™ into four parts. For every t > 0 and a £ {S', T}, 
A"(t) 

m = ~//E F b(Q*W-)/» n ) F b {QlW-)/(vn)) + F d {Ql{t?-)/n n ) - F d (Q™(t? -)./(/m)) 

^ i=l 


□ 


A"(t) 


E (W2(*"-)/M) + WS(*"-)/(A"0) - 8 

' A n (s) 


QS(*"-) 

fm 


+e Q n a (s-)d 


- / Qa(s)^ 


+ E (>(*) - /r"min(f, W"( 0 ))) . 
\/n V / 


( 8 . 112 ) 


We will show that each of these converges to zero in probability as n -> oo. 

Fix t > 0 and a € {S, T} and select arbitrary constants e,rj > 0. We first show that both 


limsupP ( -1= E \F b (Q n a (t?-)/^n)) + F b (Q n a (t2-)/^n))\ > \ ) < \ 
ra-s-oo l V n 4/4 


A n {t) 


and 


( 


A n (t ) 


limsupP -= lFd(Q2(t?-)/(»n)) + F d (Q2(t?-)/frn))l > - < {. 

oo \Vn 4 / 4 


(8.113) 


(8.114) 


We will prove (8.113) and the proof of (8.114) follows trivially. The right hand side of (8.113) can be 
expanded: 


A n (t) 

P ( Vh E \F b (QZ(t?-)/(»n))+FM(t?-)/(»n))\ > 6 
2=1 


< P A n (t) sup \F b {Ql{s)/^) + F b {Q n a {s)/^n))\ > 


(8.115) 


< 


s6[0,t] 

sup Qa(s) > Ky/n] + P(A n (t) > 2 fmt) + P ( sup I F b (sy/n / (nn)) + F b {sy/n/{fm)) I > 

, S G[0,t] / \s<K 


8fity/n J 
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Fix K > 0 so that by (5.57) of Lemma 5.1, 


sup Q”(s) > Ky/n) < 
»e[o,t] J s 


(8.116) 


Notice that for any 8 > 0 we have that 


sup 

s<K 


Sy/n 

Sy/n 

< 

Ky/n 

Ky/h 

fin 

fin 


fl n 

fin 


< 


for sufficiently large n. Set 5 = e/(16 fitd). The first result, (8.113), follows from (8.115), (8.116), and Lemmas 
3.3 and 5.3. Likewise, (8.114) follows from an identical argument. 

Next we show that 


lim sup P 

n—>oo 


i 


A n (t) 

s 


F b (Q n a (t?-)/(Lm)) + F d {Ql{t™-)/{nn)) 9 


Qm~) 


fin 




(8.117) 


We explore the derivative of the abandonment and balking distributions: 

A” 00 


— y 

- frt 1 -—^ 


i= 1 


WS(*"-)/(/"0) + F d (Q n a (t?-)/(fm)) - 9 




< 


A n (t) sup 
i se[o,t] 


F b (Q2(s)/( M n)) + F d {Ql{s)/^n)) - 9 


fin 

«(»> 


>i 




> 


ey/n 


(8.118) 


< 


sup Qa(s) > Kyfn +P (A n (t) > 2/j.nt ) 
se[o,t] J 


+1 sup 

\se[o,K/ij.] 


9s 

F b (s/y/n ) + F d {s/y/n) - -j= 

yjn 


> 


8 [ityfn 


Recall that the derivatives at zero of both F b and F d exist and sum to 9. Hence, for a any given 6 > 0 there 
exists an ho such that 

Fb{s) , F d (s) 


sup 

s<h 


<8 


for all h < ho- Let 8 = and let ho be a constant so that the above inequality holds. Now choose an no 
such that K/(fiy/n) < ho for all n > no- It follows that for each n > no, 


sup | F b (s) + F d (s ) - 9s | = sup 

s£[0,K/(iiy/n)] sE[0,K/fi] 


F{s/y/h) - 


9s 


<8 


K _ e 

Hy/n 8 [itsfn 


(8.119) 


The result (8.117) follows from (8.116), (8.118), (8.119) and Lemma 3.3. 
Third, we show that 


lim sup P I 

n—too \ 


Q n a (s-)d 


A n (s ) 
h 


- / Q n a (s)ds 


> e <T]. 


( 8 . 120 ) 


By Proposition 5.15, the process {Q|)(s),s < t} is tight. Consider a subsequence {n'} over which the process 
Qq has a limit, say Q a . By the Skorohod Representation Theorem, there exists an alternative probability 

space on which are defined a sequence {(Q& , A™ ), n > 1} and, by Lemma 3.6, a limit process ( Q a , A a ) such 
that 


(QZ,aZ)2(qz,aZ) 
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for each n' and such that (Q” ,A™ ) —► (Q a ,A a ) almost surely as n' —> oo. It is also true that 


Ql ( s-)d 


A n \s) 


D 


= / Q n As-)d 


A n \s) 


and 


Qa(s)ds= I QZ{s)ds 


for each n '. Applying Lemma 8.3 from [? ] twice we have 


sup 

uG[ 0 ,t] 


( An 1 ' ( \ \ ru . 

J - Q(s)ds 


and 


sup 

ue[o,t] 


/ Q™ ( s)ds - / Q(s)ds 

/o Jo 


almost surely as n' —> oo, so that 


sup 

nG[0,£] 


Ql («-)rf 


A n '{s) 


- Qa ( s)ds 


0 


almost surely as n' —> oo. It follows that in our original probability space 


sup 

u(E[0,£] 



-> o 


as n! —> oo. This limit holds on the arbitrarily chosen subsequence {n 7 }. Hence (8.120) holds. 

Consider the fourth term on the right hand side of (8.112). The service times of initial jobs that are 
actually served constitute an i.i.d sequence. This sequence obeys a weak law of large numbers, as does the 
renewal process, S n , constructed by these service times over intervals of time that are of the order 1 f\Jn. 
Hence, by Proposition 5.15 


sup 

S n {t) - Unt 

> ey/n ] < P ( sup 

S U (t) - H n t 

y t<w"(o) 


J \t<K/y/n 



> e 


Vn + P (W n (0) > K/y/n) < rj, (8.121) 


holds for sufficiently large n. 

Finally, our result follows from (8.112), (8.113), (8.114), (8.117), (8.120), and Lemma 3.3. O. 


8.2 Numerical Tables 
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Table 1: Simulated Results vs. Heavy Traffic Approximations Arrivals = Exp(/i + yjji ■ f3), Service 

Exp(/x), Balking = Exp(0h), Reneging = Exp(0 r ) 


Q a 

P 

P 

P 

o b 

e r 

Q 

Qrou 

W 

Wrou 

R 

Rr.OU 

B 


Brou 

Qs 

1.1 

100 

1 

.1 

,i 

49.48 ± .199 

50.74 

.481 ± .002 

.5074 

.044 ± .0002 

.0461 

.048 ± .0002 

.0461 

Qt 

1.1 

100 

1 

.1 

,i 

50.18 ± .203 

50.74 

.475 ± .002 

.5074 

.043 ± .0002 

.0461 

.049 ± .0002 

.0461 

Qs 

1.1 

25 

.5 

.1 

,i 

15.25 ± .071 

15.24 

.587 ± .003 

.6099 

.0519 ± .0002 

.0554 

.0585 ± 

.0003 

.0554 

Qt 

1.1 

25 

.5 

.1 

,i 

15.59 ± .073 

15.24 

.578 ± .003 

.6099 

.0510 ± .0002 

.0554 

.0597 ± 

.0003 

.0554 

Qs 

1.1 

4 

.2 

.1 

,i 

4.47 ± .021 

4.40 

1.04 ± .0006 

1.100 

.0813 ± .0004 

.1000 

.102 ± .0006 

.1000 

Qt 

1.1 

4 

.2 

.1 

,i 

4.67 ± .023 

4.40 

1.02 ± .0006 

1.100 

.0791 ± .0004 

.1000 

.106 ± .0006 

.1000 

Qs 

1.01 

100 

.1 

.1 

,i 

19.79 ± .094 

19.78 

.195 ± .001 

.1978 

.0186 ± .0001 

.0195 

.0195 ± 

.0001 

.0195 

Qt 

1.01 

100 

.1 

.1 

,i 

19.98 ± .095 

19.78 

.194 ± .001 

.1978 

.0186 ± .0001 

.0195 

.0196 ± 

.0001 

.0195 

Qs 

1.01 

25 

.05 

.1 

,i 

9.34 ± .045 

9.39 

.363 ± .0018 

.3956 

.0333 ± .0002 

.03718 

.0362 ± 

.0002 

.03718 

Qt 

1.01 

25 

.05 

.1 

,i 

9.51 ± .046 

9.39 

.360 ± .0018 

.3956 

.0330 ± .0002 

.03718 

.0368 ± 

.0002 

.03718 

Qs 

1.01 

4 

.02 

.1 

,i 

3.61 ± .019 

3.64 

.855 ± .005 

.9104 

.0681 ± .0005 

.0901 

.0843 ± 

.0006 

.0901 

Qt 

1.01 

4 

.02 

.1 

,i 

3.77 ± .020 

3.64 

.842 ± .005 

.9104 

.0667 ± .0004 

.0901 

.0873 ± 

.0006 

.0901 

Qs 

1 

100 

0 

.1 

,i 

17.82 ± .083 

17.84 

.176 ± .0009 

.1784 

.0168 ± .0001 

.0178 

.0176 ± 

.0001 

.0178 

Qt 

1 

100 

0 

.1 

,i 

17.98 ± .085 

17.84 

.175 ± .0009 

.1784 

.0168 ± .0001 

.0178 

.0177 ± 

.0001 

.0178 

Qs 

1 

25 

0 

.1 

,i 

8.86 ± .046 

8.92 

.345 ± .002 

.3568 

.0317 ± .0002 

.0356 

.0344 ± 

.0002 

.0356 

Qt 

1 

25 

0 

.1 

,i 

9.03 ± .047 

8.92 

.343 ± .002 

.3568 

.0314 ± .0002 

.0356 

.0350 ± 

.0002 

.0356 

Qs 

1 

4 

0 

.1 

,i 

3.52 ± .019 

3.56 

.832 ± .005 

.8920 

.0670 ± .0005 

.0892 

.0817 ± 

.0005 

.0892 

Qt 

1 

4 

0 

.1 

,i 

3.67 ± .021 

3.56 

.820 ± .005 

.8920 

.0657 ± .0005 

.0892 

.0848 ± 

.0005 

.0892 

Qs 

.99 

100 

-.1 

.1 

,i 

15.95 ± .067 

16.14 

.157 ± .0007 

.1614 

.0152 ± .0001 

.0163 

.0157 ± 

.0001 

.0163 

Qt 

.99 

100 

-.1 

.1 

,i 

16.10 ± .068 

16.14 

.157 ± .0007 

.1614 

.0151 ± .0001 

.0163 

.0159 ± 

.0001 

.0163 

Qs 

.99 

25 

-.05 

.1 

,i 

8.40 ± .046 

8.48 

.328 ± .002 

.3392 

.0302 ± .0002 

.0342 

.0326 ± 

.0002 

.0342 

Qt 

.99 

25 

-.05 

.1 

,i 

8.55 ± .046 

8.48 

.325 ± .002 

.3392 

.0299 ± .0002 

.0342 

.0332 ± 

.0002 

.0342 

Qs 

.99 

4 

-.02 

.1 

,i 

3.44 ± .020 

3.49 

.814 ± .0058 

.8741 

.0657 ± .0005 

.0882 

.0801 ± 

.0005 

.0882 

Qt 

.99 

4 

-.02 

.1 

,i 

3.58 ± .022 

3.49 

.803 ± .0055 

.8741 

.0644 ± .0005 

.0882 

.0830 ± 

.0005 

.0882 

Qs 

.9 

100 

-1 

.1 

,i 

7.12 ± .0264 

7.77 

.0708 ± .0003 

.0777 

.00692 ± .0001 

.0086 

.0071 ± 

.0001 

.0086 

Qt 

.9 

100 

-1 

.1 

,i 

7.16 ± .0267 

7.77 

.0707 ± .0003 

.0777 

.00691 ± .0001 

.0086 

.0071 ± 

.0001 

.0086 

Qs 

.9 

25 

-.5 

.1 

,i 

5.19 ± .024 

5.61 

.204 ± .001 

.2246 

.0192 ± .0001 

.0249 

.0204 ± 

.0001 

.0249 

Qt 

.9 

25 

-.5 

.1 

,i 

5.27 ± .024 

5.61 

.203 ± .001 

.2246 

.0192 ± .0001 

.0249 

.0207 ± 

.0001 

.0249 

Qs 

.9 

4 

-.2 

.1 

,i 

2.74 ± .013 

2.93 

.653 ± .0036 

.7328 

.0542 ± .0004 

.0814 

.0645 ± 

.0005 

.0814 

Qt 

.9 

4 

-.2 

.1 

,i 

2.84 ± .014 

2.93 

.647 ± .0035 

.7328 

.0534 ± .0004 

.0814 

.0667 ± 

.0005 

.0814 

Qs 

.8 

100 

-2 

.1 

,i 

3.73 ± .0098 

4.59 

.0372 ± .0001 

.0459 

.0037 ± .0001 

.0057 

.0037 ± 

.0001 

.0057 

Qt 

.8 

100 

-2 

.1 

,i 

3.75 ± .0099 

4.59 

.0372 ± .0001 

.0459 

.0037 ± .0001 

.0057 

.0037 ± 

.0001 

.0057 

Qs 

.8 

25 

-1 

.1 

,i 

3.22 ± .011 

3.88 

.127 ± .0005 

.1555 

.0122 ± .0001 

.0194 

.0127 ± 

.0001 

.0194 

Qt 

.8 

25 

-1 

.1 

,i 

3.25 ± .012 

3.88 

.127 ± .0005 

.1555 

.0122 ± .0001 

.0194 

.0129 ± 

.0001 

.0194 

Qs 

.8 

4 

-.4 

.1 

,i 

2.11 ± .012 

2.44 

.508 ± .0036 

.6113 

.0432 ± .0003 

.0764 

.0496 ± 

.0003 

.0764 

Qt 

.8 

4 

-.4 

.1 

,i 

2.18 ± .013 

2.44 

.505 ± .0035 

.6113 

.0428 ± .0003 

.0764 

.0511 ± 

.0003 

.0764 
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Table 2: Simulated Results vs. Heavy Traffic Approximations Arrivals = Exp(/i + v fji • /?), Service = 

LogNormal(l//x,l/V 2 )> 

Balking = Unifonn(O,l/0b), Reneging = Unifonn(0,l/#r) 


Q a 

P 

P 

P 

e b 

e r 

Q 

Qrou 

W 

Wnou 

R 

Rrou 

B 


Brou 

Qs 

1.1 

100 

1 

.1 

,i 

48.17 ± .174 

50.74 

.468 ± .002 

.5074 

.044 ± .0002 

.0461 

.048 ± .0002 

.0461 

Qt 

1.1 

100 

1 

.1 

,i 

48.85 ± .177 

50.74 

.462 ± .002 

.5074 

.043 ± .0002 

.0461 

.049 ± .0002 

.0461 

Qs 

1.1 

25 

.5 

.1 

,i 

14.65 ± .073 

15.24 

.563 ± .003 

.6099 

.0517 ± .0002 

.0554 

.0587 ± 

.0003 

.0554 

Qt 

1.1 

25 

.5 

.1 

,i 

14.97 ± .076 

15.24 

.555 ± .003 

.6099 

.0508 ± .0002 

.0554 

.0599 ± 

.0003 

.0554 

Qs 

1.1 

4 

.2 

.1 

,i 

4.17 ± .022 

4.40 

.968 ± .0006 

1.100 

.0813 ± .0004 

.1000 

.104 ± .0006 

.1000 

Qt 

1.1 

4 

.2 

.1 

,i 

4.35 ± .024 

4.40 

.946 ± .0006 

1.100 

.0786 ± .0004 

.1000 

.108 ± .0006 

.1000 

Qs 

1.01 

100 

.1 

.1 

,i 

19.35 ± .086 

19.78 

.191 ± .0009 

.1978 

.0185 ± .0001 

.0195 

.0194 ± 

.0001 

.0195 

Qt 

1.01 

100 

.1 

.1 

,i 

19.55 ± .087 

19.78 

.189 ± .0009 

.1978 

.0184 ± .0001 

.0195 

.0196 ± 

.0001 

.0195 

Qs 

1.01 

25 

.05 

.1 

,i 

9.04 ± .047 

9.39 

.352 ± .0020 

.3956 

.0331 ± .0002 

.03718 

.0362 ± 

.0002 

.03718 

Qt 

1.01 

25 

.05 

.1 

,i 

9.22 ± .048 

9.39 

.349 ± .00189 

.3956 

.0328 ± .0002 

.03718 

.0369 ± 

.0002 

.03718 

Qs 

1.01 

4 

.02 

.1 

,i 

3.40 ± .019 

3.64 

.795 ± .006 

.9104 

.0684 ± .0005 

.0901 

.0849 ± 

.0005 

.0901 

Qt 

1.01 

4 

.02 

.1 

,i 

3.54 ± .020 

3.64 

.792 ± .005 

.9104 

.0667 ± .0005 

.0901 

.0884 ± 

.0006 

.0901 

Qs 

1 

100 

0 

.1 

,i 

17.45 ± .082 

17.84 

.172 ± .0009 

.1784 

.0167 ± .0001 

.0178 

.0174 ± 

.0001 

.0178 

Qt 

1 

100 

0 

.1 

,i 

17.62 ± .084 

17.84 

.171 ± .0009 

.1784 

.0167 ± .0001 

.0178 

.0176 ± 

.0001 

.0178 

Qs 

1 

25 

0 

.1 

,i 

8.62 ± .037 

8.92 

.335 ± .002 

.3568 

.0317 ± .0002 

.0356 

.0344 ± 

.0002 

.0356 

Qt 

1 

25 

0 

.1 

,i 

8.78 ± .038 

8.92 

.333 ± .002 

.3568 

.0314 ± .0002 

.0356 

.0351 ± 

.0002 

.0356 

Qs 

1 

4 

0 

.1 

,i 

3.30 ± .016 

3.56 

.773 ± .006 

.8920 

.0663 ± .0005 

.0892 

.0825 ± 

.0005 

.0892 

Qt 

1 

4 

0 

.1 

,i 

3.44 ± .017 

3.56 

.761 ± .005 

.8920 

.0647 ± .0005 

.0892 

.0860 ± 

.0005 

.0892 

Qs 

.99 

100 

-.1 

.1 

,i 

15.78 ± .077 

16.14 

.156 ± .0008 

.1614 

.0152 ± .0001 

.0163 

.0157 ± 

.0001 

.0163 

Qt 

.99 

100 

-.1 

.1 

,i 

15.92 ± .078 

16.14 

.155 ± .0008 

.1614 

.0151 ± .0001 

.0163 

.0159 ± 

.0001 

.0163 

Qs 

.99 

25 

-.05 

.1 

,i 

8.18 ± .040 

8.48 

.319 ± .002 

.3392 

.0301 ± .0002 

.0342 

.0327 ± 

.0002 

.0342 

Qt 

.99 

25 

-.05 

.1 

,i 

8.32 ± .042 

8.48 

.317 ± .002 

.3392 

.0298 ± .0002 

.0342 

.0333 ± 

.0002 

.0342 

Qs 

.99 

4 

-.02 

.1 

,i 

3.24 ± .018 

3.49 

.761 ± .006 

.8741 

.0655 ± .0005 

.0882 

.0812 ± 

.0005 

.0882 

Qt 

.99 

4 

-.02 

.1 

,i 

3.38 ± .019 

3.49 

.750 ± .006 

.8741 

.0640 ± .0004 

.0882 

.0846 ± 

.0005 

.0882 

Qs 

.9 

100 

-1 

.1 

,i 

7.09 ± .0252 

7.77 

.0705 ± .0003 

.0777 

.00694 ± .0001 

.0086 

.0071 ± 

.0001 

.0086 

Qt 

.9 

100 

-1 

.1 

,i 

7.13 ± .0256 

7.77 

.0704 ± .0003 

.0777 

.00693 ± .0001 

.0086 

.0071 ± 

.0001 

.0086 

Qs 

.9 

25 

-.5 

.1 

,i 

5.10 ± .027 

5.61 

.200 ± .001 

.2246 

.0193 ± .0001 

.0249 

.0205 ± 

.0001 

.0249 

Qt 

.9 

25 

-.5 

.1 

,i 

5.18 ± .028 

5.61 

.200 ± .001 

.2246 

.0192 ± .0001 

.0249 

.0207 ± 

.0001 

.0249 

Qs 

.9 

4 

-.2 

.1 

,i 

2.61 ± .015 

2.93 

.619 ± .004 

.7328 

.0541 ± .0004 

.0814 

.0655 ± 

.0005 

.0814 

Qt 

.9 

4 

-.2 

.1 

,i 

2.71 ± .016 

2.93 

.612 ± .004 

.7328 

.0532 ± .0004 

.0814 

.0680 ± 

.0005 

.0814 

Qs 

.8 

100 

-2 

.1 

,i 

3.72 ± .0098 

4.59 

.0371 ± .0001 

.0459 

.0037 ± .0001 

.0057 

.0037 ± 

.0001 

.0057 

Qt 

.8 

100 

-2 

.1 

,i 

3.74 ± .0099 

4.59 

.0371 ± .0001 

.0459 

.0037 ± .0001 

.0057 

.0037 ± 

.0001 

.0057 

Qs 

.8 

25 

-1 

.1 

,i 

3.19 ± .013 

3.88 

.126 ± .0006 

.1555 

.0123 ± .0001 

.0194 

.0127 ± 

.0001 

.0194 

Qt 

.8 

25 

-1 

.1 

,i 

3.22 ± .014 

3.88 

.126 ± .0006 

.1555 

.0122 ± .0001 

.0194 

.0129 ± 

.0001 

.0194 

Qs 

.8 

4 

-.4 

.1 

,i 

2.02 ± .011 

2.44 

.481 ± .0034 

.6113 

.0429 ± .0004 

.0764 

.0505 ± 

.0004 

.0764 

Qt 

.8 

4 

-.4 

.1 

,i 

2.08 ± .011 

2.44 

.478 ± .0033 

.6113 

.0424 ± .0004 

.0764 

.0522 ± 

.0004 

.0764 
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Table 3: Simulated Results vs. Heavy Traffic Approximations Arrivals = Exp(/i + v /Jl ■ /?), Service 

Exp(/x), Balking = Exp(0h), Reneging = Exp(0 r ) 


Q a 

P 

P 

P 

o b 

e r 

Q 

Qrou 

W 

WroU 

R 

Rr.OU 

B 

Brou 

Qs 

1.1 

100 

1 

1 

i 

7.96 ± .012 

7.88 

.075 ± 0001 

.0788 

.0637 ± .0001 

.0717 

.075 ± .0001 

.0717 

Qt 

1.1 

100 

1 

1 

i 

8.21 ± .013 

7.88 

.074 ± .0001 

.0788 

.0624 ± .0001 

.0717 

.077 it .0001 

.0717 

Qs 

1.1 

25 

.5 

1 

i 

3.38 it : .006 

3.32 

.125 ± .0002 

.133 

.122 ± .0001 

.1209 

.091 ± .0002 

.1209 

Qt 

1.1 

25 

.5 

1 

i 

3.56 ± .006 

3.32 

.122 ± .0003 

.133 

.0088 ± .0002 

.1209 

.127 st .0002 

.1209 

Qs 

1.1 

4 

.2 

1 

i 

1.24 ± .003 

1.20 

.277 ± .0009 

.301 

.124 ± .0004 

.2736 

.240 ± .0005 

.2736 

Qt 

1.1 

4 

.2 

1 

i 

1.37 ± .003 

1.20 

.268 ± .0009 

.301 

.117 ± .0004 

.2736 

.255 ± .0005 

.2736 

Qs 

1.01 

100 

.1 

1 

i 

5.78 ± .009 

5.82 

.055 =t .0001 

.0582 

.0483 ± .0001 

.0576 

.0552 ± .0001 

.0576 

Qt 

1.01 

100 

.1 

1 

i 

5.94 ± .009 

5.82 

.055 =t .0001 

.0582 

.0476 ± .0001 

.0576 

.0567 ± .0001 

.0576 

Qs 

1.01 

25 

.05 

1 

i 

2.83 ± .005 

2.86 

.105 ± .0002 

.1146 

.079 it .0002 

.1135 

.103 it .0002 

.1135 

Qt 

1.01 

25 

.05 

1 

i 

2.97 ± .005 

2.86 

.103 ± .0002 

.1146 

.077 ± .0002 

.1135 

.107 it .0002 

.1135 

Qs 

1.01 

4 

.02 

1 

i 

1.14 ± .003 

1.13 

.255 ± .0009 

.2839 

.118 it .0005 

.2811 

.222 ± .0006 

.2811 

Qt 

1.01 

4 

.02 

1 

i 

1.25 ± .003 

1.13 

.249 ± .0009 

.2839 

.112 ± .0004 

.2811 

.235 ± .0006 

.2811 

Qs 

1 

100 

0 

1 

i 

5.58 ± .009 

5.64 

.0537 ± .0001 

.0564 

.047 ± .0001 

.0564 

.053 ± .0001 

.0564 

Qt 

1 

100 

0 

1 

i 

5.74 ± .009 

5.64 

.0531 ± .0001 

.0564 

.046 ± .0001 

.0564 

.055 ± .0001 

.0564 

Qs 

1 

25 

0 

1 

i 

2.78 ± .005 

2.82 

.104 ± .0002 

.1128 

.078 ± .0002 

.1128 

.102 it .0002 

.1128 

Qt 

1 

25 

0 

1 

i 

2.92 ± .005 

2.82 

.102 ± .0002 

.1128 

.076 ± .0002 

.1128 

.106 it .0002 

.1128 

Qs 

1 

4 

0 

1 

i 

1.12 it .003 

1.12 

.253 it .0010 

.2820 

.117 ± .0005 

.2820 

.220 ± .0006 

.2820 

Qt 

1 

4 

0 

1 

i 

1.23 ± .003 

1.12 

.247 ± .0009 

.2820 

.110 ± .0004 

.2820 

.233 ± .0006 

.2820 

Qs 

.99 

100 

-.1 

1 

i 

5.39 ± .009 

5.46 

.052 ± .0001 

.0551 

.045 ± .0001 

.0546 

.052 ± .0001 

.0546 

Qt 

.99 

100 

-.1 

1 

i 

5.54 ± .009 

5.46 

.051 ± .0001 

.0551 

.045 it .0001 

.0546 

.053 ± .0001 

.0546 

Qs 

.99 

25 

-.05 

1 

i 

2.73 ± .005 

2.77 

.102 ± .0002 

.1110 

.077 ± .0002 

.1121 

.100 ± .0002 

.1121 

Qt 

.99 

25 

-.05 

1 

i 

2.86 ± .006 

2.77 

.100 ± .0002 

.1110 

.075 ± .0002 

.1121 

.104 ± .0002 

.1121 

Qs 

.99 

4 

-.02 

1 

i 

1.11 ± .003 

1.12 

.251 it : .0009 

.2802 

.117 ± .0005 

.2831 

.219 ± .0006 

.2831 

Qt 

.99 

4 

-.02 

1 

i 

1.22 ± .003 

1.12 

.245 ± .0009 

.2802 

.110 st .0005 

.2831 

.232 ± .0006 

.2831 

Qs 

.9 

100 

-1 

1 

i 

3.88 ± .005 

4.16 

.0377 ± .0001 

.0416 

.034 ± .0001 

.0462 

.037 ± .0001 

.0462 

Qt 

.9 

100 

-1 

1 

i 

3.97 ± .006 

4.16 

.0375 ± .0001 

.0416 

.034 * .0001 

.0462 

.038 ± .0001 

.0462 

Qs 

.9 

25 

-.5 

1 

i 

2.25 ± .005 

2.41 

.085 ± .0002 

.0964 

.066 ± .0002 

.1071 

.083 ± .0002 

.1071 

Qt 

.9 

25 

-.5 

1 

i 

2.36 ± .005 

2.41 

084 ± .0002 

.0964 

.065 it .0002 

.1071 

.087 ± .0002 

.1071 

Qs 

.9 

4 

-.2 

1 

i 

1.01 ± .003 

1.05 

.229 ± .0009 

.2646 

.109 it .0005 

.2940 

.200 ± .0006 

.2940 

Qt 

.9 

4 

-.2 

1 

i 

1.10 ± .003 

1.05 

.224 ± .0009 

.2646 

.103 it .0004 

.2940 

.211 ± .0007 

.2940 

Qs 

.8 

100 

-2 

1 

i 

2.70 ± .005 

3.19 

.026 ± .0001 

.0319 

.024 ± .0001 

.0399 

.026 ± .0001 

.0399 

Qt 

.8 

100 

-2 

1 

i 

2.75 ± .005 

3.19 

.026 ± .0001 

.0319 

.024 ± .0001 

.0399 

.027 ± .0001 

.0399 

Qs 

.8 

25 

-1 

1 

i 

1.80 =t .003 

2.08 

.069 ± .0002 

.0832 

.055 ± .0002 

.1040 

.067 it .0002 

.1040 

Qt 

.8 

25 

-1 

1 

i 

1.87 ± .003 

2.08 

.068 it .0002 

.0832 

.054 ± .0002 

.1040 

.069 ± .0002 

.1040 

Qs 

.8 

4 

-.4 

1 

i 

.89 ± .002 

.9947 

.204 ± .0008 

.2486 

.100 ± .0005 

.3108 

.178 ± .0007 

.3108 

Qt 

.8 

4 

-.4 

1 

i 

.97 ± .002 

.9947 

.200 ± .0008 

.2486 

.095 ± .0004 

.3108 

.188 ± .0007 

.3108 
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Table 4: Simulated Results vs. Heavy Traffic Approximations Arrivals = Exp(/i + v fji • /?), Service = 

LogNormal(l//x,l/V 2 )> 

Balking = Unifonn(O,l/0b), Reneging = Unifonn(0,l/#r) 


Q a 

P 

P 

P 

e b 

e r 

Q 

Qrou 

W 

Wrou 

R 

Rrou 

B 

Rrou 

Qs 

1.1 

100 

1 

i 

i 

7.56 ± .011 

7.88 

.072 ± 0001 

.0788 

.0636 ± .0001 

.0717 

.076 ± .0001 

.0717 

Qt 

1.1 

100 

1 

i 

i 

7.79 ± .012 

7.88 

.070 ± .0001 

.0788 

.0622 ± .0001 

.0717 

.078 ± .0001 

.0717 

Qs 

1.1 

25 

.5 

i 

i 

3.11 ± .005 

3.32 

.114 ± .0003 

.133 

.091 ± .0001 

.1209 

.124 ± .0002 

.1209 

Qt 

1.1 

25 

.5 

i 

i 

3.27 ± .006 

3.32 

.111 ± .0002 

.133 

.088 ± .0002 

.1209 

.131 ± .0002 

.1209 

Qs 

1.1 

4 

.2 

i 

i 

1.06 ± .002 

1.20 

.232 ± .001 

.301 

.118 ± .0005 

.2736 

.266 ± .0006 

.2736 

Qt 

1.1 

4 

.2 

i 

i 

1.15 ± .002 

1.20 

.225 ± .001 

.301 

.105 ± .0004 

.2736 

.289 ± .0007 

.2736 

Qs 

1.01 

100 

.1 

i 

i 

5.57 ± .009 

5.82 

.052 ± .0001 

.0582 

.0470 ± .0001 

.0576 

.0557 ± .0001 

.0576 

Qt 

1.01 

100 

.1 

i 

i 

5.73 ± .009 

5.82 

.051 ± .0001 

.0582 

.0463 ± .0001 

.0576 

.0573 ± .0001 

.0576 

Qs 

1.01 

25 

.05 

i 

i 

2.65 ± .004 

2.86 

.098 ± .0002 

.1146 

.080 ± .0002 

.1135 

.106 ± .0002 

.1135 

Qt 

1.01 

25 

.05 

i 

i 

2.78 ± .005 

2.86 

.096 ± .0002 

.1146 

.078 ± .0002 

.1135 

.111 ± .0002 

.1135 

Qs 

1.01 

4 

.02 

i 

i 

.99 ± .002 

1.13 

.218 ± .001 

.2839 

.113 ± .0005 

.2811 

.247 ± .0006 

.2811 

Qt 

1.01 

4 

.02 

i 

i 

1.07 ± .002 

1.13 

.211 ± .001 

.2839 

.101 ± .0004 

.2811 

.268 ± .0007 

.2811 

Qs 

1 

100 

0 

i 

i 

5.38 ± .008 

5.64 

.0516 ± .0001 

.0564 

.047 ± .0001 

.0564 

.054 ± .0001 

.0564 

Qt 

1 

100 

0 

i 

i 

5.53 ± .009 

5.64 

.0510 ± .0001 

.0564 

.046 ± .0001 

.0564 

.055 ± .0001 

.0564 

Qs 

1 

25 

0 

i 

i 

2.59 ± .005 

2.82 

.095 ± .0003 

.1128 

.079 ± .0002 

.1128 

.104 ± .0002 

.1128 

Qt 

1 

25 

0 

i 

i 

2.72 ± .005 

2.82 

.094 ± .0002 

.1128 

.077 ± .0002 

.1128 

.109 ± .0002 

.1128 

Qs 

1 

4 

0 

i 

i 

.98 ± .002 

1.12 

.217 =t .001 

.2820 

.113 ± .0005 

.2820 

.245 ± .0006 

.2820 

Qt 

1 

4 

0 

i 

i 

1.06 ± .002 

1.12 

.210 ± .001 

.2820 

.100 ± .0004 

.2820 

.265 ± .0006 

.2820 

Qs 

.99 

100 

-.1 

i 

i 

5.18 ± .008 

5.46 

.049 ± .0001 

.0551 

.045 ± .0001 

.0546 

.052 ± .0001 

.0546 

Qt 

.99 

100 

-.1 

i 

i 

5.33 ± .008 

5.46 

.049 ± .0001 

.0551 

.045 ± .0001 

.0546 

.053 ± .0001 

.0546 

Qs 

.99 

25 

-.05 

i 

i 

2.54 ± .005 

2.77 

.094 ± .0002 

.1110 

.078 ± .0002 

.1121 

.102 ± .0002 

.1121 

Qt 

.99 

25 

-.05 

i 

i 

2.67 ± .005 

2.77 

.093 ± .0002 

.1110 

.075 ± .0002 

.1121 

.108 ± .0002 

.1121 

Qs 

.99 

4 

-.02 

i 

i 

.97 ± .002 

1.12 

.214 ± .0009 

.2802 

.112 ± .0005 

.2831 

.242 ± .0006 

.2831 

Qt 

.99 

4 

-.02 

i 

i 

1.05 ± .002 

1.12 

.208 ± .0009 

.2802 

.101 ± .0005 

.2831 

.263 ± .0006 

.2831 

Qs 

.9 

100 

-1 

i 

i 

3.77 ± .007 

4.16 

.0364 ± .0001 

.0416 

.034 ± .0001 

.0462 

.038 ± .0001 

.0462 

Qt 

.9 

100 

-1 

i 

i 

3.86 ± .007 

4.16 

.0346 ± .0001 

.0416 

.034 ± .0001 

.0462 

.039 ± .0001 

.0462 

Qs 

.9 

25 

-.5 

i 

i 

2.13 ± .004 

2.41 

.080 ± .0002 

.0964 

.068 ± .0002 

.1071 

.085 ± .0002 

.1071 

Qt 

.9 

25 

-.5 

i 

i 

2.23 ± .004 

2.41 

.078 ± .0002 

.0964 

.065 ± .0002 

.1071 

.089 ± .0002 

.1071 

Qs 

.9 

4 

-.2 

i 

i 

.89 ± .002 

1.05 

.198 ± .0009 

.2646 

.106 ± .0005 

.2940 

.223 ± .0006 

.2940 

Qt 

.9 

4 

-.2 

i 

i 

.96 ± .002 

1.05 

.192 ± .0009 

.2646 

.096 ± .0004 

.2940 

.241 ± .0007 

.2940 

Qs 

.8 

100 

-2 

i 

i 

2.65 ± .004 

3.19 

.026 ± .0001 

.0319 

.024 ± .0001 

.0399 

.026 ± .0001 

.0399 

Qt 

.8 

100 

-2 

i 

i 

2.70 ± .004 

3.19 

.026 ± .0001 

.0319 

.024 ± .0001 

.0399 

.027 ± .0001 

.0399 

Qs 

.8 

25 

-1 

i 

i 

1.72 ± .003 

2.08 

.065 ± .0002 

.0832 

.056 ± .0002 

.1040 

.069 ± .0002 

.1040 

Qt 

.8 

25 

-1 

i 

i 

1.79 ± .003 

2.08 

.064 ± .0002 

.0832 

.055 ± .0002 

.1040 

.072 ± .0002 

.1040 

Qs 

.8 

4 

-.4 

i 

i 

.79 ± .002 

.9947 

.180 ± .0008 

.2486 

.090 ± .0005 

.3108 

.199 ± .0007 

.3108 

Qt 

.8 

4 

-.4 

i 

i 

.86 ± .002 

.9947 

.176 ± .0008 

.2486 

.099 ± .0004 

.3108 

.215 ± .0007 

.3108 
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